AD-A188  2M  C 


1/1 


UNCLASSIFIED 


r/c  / 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  Of  STANDARDS  196/  A 


.v^ftON  OF  THIS  PAGE 


14  REPORT  SECURITY  CLASSIFICATION 


Unclassified 


24  SECURITY  CLASSIFICATION  AUTHORU 


2b  DECLASSIFICATION  /  OOWNG 


4  PERFORMING  ORGANIZATION  R 


AD- A 188  204 


4RKINGS 


iVAILABIUTY  OF  REPORT 
Approved  for  public  release 
distribution  unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUNMEBairy 

AFOSR-TR-  87-10  57 


6*.  NAME  OF  PERFORMING  ORGANIZATION 

Northeastern  University 

sMffice  symbol 

(If  epplkeb/e) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research/NL 

6c.  ADDRESS  (Cfy.  Stote.  end  ZIP  Code) 

Office  of  Sponsored  Research 

Northeastern  University 

_ Boston.  MA  02115 _ 

7b  AOORESS  (City.  State,  end  ZIP  Code) 

Building  410 

Bolling  AFB,  D.C.  20332-6448 

8*.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

AFOSR 

8b  OFFICE  SYMBOL 
(If  oppliceble) 

NL 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F49620-87-C-0013 

Be  ADDRESS  (City.  Stott,  end  ZIP  Code) 
Building  410 
Bolling  AFB 

2-644 


11  TITLE  (Include  Security  Geaificetion ) 

Center  for  the  Study  of  Rhythmic  Processes 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


MS! 


12  PERSONAL  AUTHOR(S) 
N.  Kopell,  A.H. 


134.  TYPE  OF  REPORT 

Annual 


Cohen,  E.  Marder,  K.  Sigvardt 


13b  TIME  COVEREO  114  DATE  OF  REPOI 

FROM  86/10/1  TO 87/ 10/ 1|  1987/10/20 


14  DATE  OF  REPORT  (Year,  Mont*.  Day) 

1987/10/20 


COSATI  COOES 


GROUP  SUB-GROUP 


U  SUBJECT  TERMS  (Confirm*  on  reverie  if  necessory  end  identify  by  block  number) 

Oscillators  Central  pattern  generators  Neural  network 

Spinal  cord  Mathematical  modeling  Neuromodulators 

Regeneration  Sensory  feedback 


19  ABSTRACT  ( Continue  on  reverse  if  necetury  end  identify  by  block  number) 

The  Center  for  the  Study  of  Rhythmic  Processes  began  operation  in  the  academic  year 
1986-87.  It  brought  together  mathematicians  and  biologists  to  work  on  problems  of  mutual 
interest.  During  the  first  year,  the  Center  focused  on  two  sets  of  problems.  The  first 
involved  the  structure  and  function  of  the  intersegmental  coordinating  system  of  the  verte¬ 
brate  spinal  cord,  for  which  the  lamprey  provides  an  excellent  model  system.  A  broadly 
applicable  mathematical  framework  was  constructed  for  analyzing  this  system,  and  new  mathe¬ 
matical  techniques  were  invented.  A  new  technology  was  put  into  use  which  has  the  potential 
of  providing  information  not  previously  obtainable.  Many  of  the  major  researchers  working  on 
this  preparation  were  consolidated  under  the  auspices  of  the,  center,  and  many  new  collabora¬ 
tions  were  formed.  The  other  involved  small  neural  networks,  such  as  the  stomatogastric 
ganglion  of  the  lobster.  Detailed  mechanisms  of  function  and  control  were  investigated, 
and  mathematical  tools  were  applied  to  investigate  how  the  circuits  change  under  modulation. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

DP  UNCLASSIFIED/UNLIMITED  E  SAME  AS  RPT  Q  OTIC  USERS 


22*  NAME  OF  RESPONSIBLE  INDIVIDUAL 
.  Alfred  Fregly,  Ph.D. 


00  FORM  1 473, 84  MAR  83  APR  edition  may  be  used  until  eihauited 

All  other  editions  ere  obsolete. 


21  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TELEPHONE  (Include  Aree  Code ) 
202-767-5021 


SECURITY  CLASSIFICATION  OF  ’HIS  PAGE 
UNCLASSIFIED 


AFOSR.TR.  8  7  -  1  657 


Center  for  the  Study  of  Rhythmic  Processes 
Progress  Report  -  Year  1  (10/1/86  -  10/1/87) 


Contents 


I.  Personnel  of  the  Center 

II.  Stateaent  of  Work 

III.  Mathematics  accomplished  by  Ermentrout  and  Kopell  (by  N.  Kopell) 

IV.  Work  accomplished  in  the  Marder  Laboratory  (by  E.  Marder) 

V.  Work  accomplished  in  the  Cohen  Laboratory  (by  A.  Cohen) 

VI.  Work  accomplished  by  Sigvardt  and  Williams  (by  K.  Sigvardt) 

VII.  Coordinating  activities  of  the  Center  (by  N.  Kopell) 

VIII.  Appendix  -  Workshop  Agenda 


Personnel  of  the  Center 


1.  Nancy  Kopell,  Department  of  Mathematics,  Boston  University,  P.I. 

Steven  Strogatz,  Postdoctoral  Fellow 
Stephane  Laederich,  Graduate  Student 

2.  G.  Bard  Brmentrout,  Department  of  Mathematics,  University  of  Pittsburgh 

3.  Eve  Marder,  Department  of  Biology,  Brandeis  University 

Pierre  Meyrand,  Postdoctoral  Fellow 
Robert  Zarum,  Undergraduate 
James  Weimann,  Graduate  Student 

4.  Avis  Cohen,  Neurobiology  and  Behavior,  Cornell  University 

Nicholis  Mellen,  Graduate  Student 

Tamara  Dobrof,  Graduate  Student 

Margaret  Baker,  Technician  and  Collaborator 

5.  Karen  Sigvardt,  Department  of  Neurology,  University  of  California,  Davis 

6.  Thelma  Williams,  St.  Georges  Hospital  Medical  School  (London) 

Simon  Alford,  Graduate  Student  (essentially  finished) 

7.  Andrew  McClellan,  Department  of  Physiology  and  Biophysics,  University  of 

Iowa 

8.  Irving  Epstein,  Department  of  Chemistry,  Brandeis  University 


1.  Carl  Rovainen,  Department  of  Physiology  and  Biophysics,  Washington 

University  Medical  School 

2.  Simon  Giszter,  Bizzi  Laboratory,  Brain  and  Cognitive  Sciences,  MIT 


CpBSMUWlta 

1.  Lawrence  Cohen,  Dept,  of  Physiology,  Yale  University 

2.  Humberto  Maturana,  Department  of  Biology,  University  of  Chile 

3.  Peter  Brodfeuhrer,  Neurobiology  and  Behavior,  Cornell  University 


II 


During  the  past  year,  we  have  established  a  Center  for  the  Study  of 
Rhythsic  Processes.  The  Center  brings  together  aatheaat iciana  and  biologists 
in  seainars,  workshops  and  an  ataosphere  of  intense  discussion  to  create  an 
environaent  in  which  p rob leas  of  autual  interest  can  be  aost  effectively  and 
creatively  investigated.  The  aatheaatical  work  entails  the  use  of  dynaaical 
systeas  theory  to  develop  a  broadly  applicable  theoretical  fraaework  for 
explaining  the  behavior  of  single  oscillators  as  well  as  groups  of  coupled 
oscillators.  The  physiology  explores  the  structure  of  neural  networks  that 
can  generate  oscillatory  outputs.  Insights  froa  the  aat hematics  are  used  to 
help  foraulate  productive  experiaents,  as  well  as  to  aid  in  the  interpretation 
of  the  results.  The  results  of  the  experiaents  are  used  to  infora  the 
development  of  the  aatheaatical  fraaework. 

During  the  first  year,  the  Center  focussed  on  two  sets  of  problems.  One 
involves  the  structure  and  function  of  the  intersegaental  coordinating  systea 
of  the  vertebrate  spinal  cord;  the  lamprey  provides  an  excellent  model  systea 
for  such  an  investigation.  New  aatheaatical  techniques  have  been  invented, 
and  new  technology  has  been  put  into  use  which  has  the  potential  of  providing 
information  not  previouly  obtainable.  Many  of  the  aajor  researchers  working 
on  this  preparation  have  been  consolidated  under  the  auspices  of  the  Center, 
and  aany  new  collaborations  have  formed. 

The  other  problem  involves  small  neural  networks,  such  as  are  found  in 
some  invertebrate  preparations.  Detailed  mechanisms  of  function  and  control 
of  neural  oscillators  are  being  investigated  in  a  system  in  which  there  are 
few  cells  and  the  circuitry  is  near  to  being  understood.  Mathematical  tools 
are  being  applied  to  investigate  how  the  circuits  change  under  modulation. 
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HI.  Mathematics  done  by  Braentrout  and  Kooell 

Much  progress  has  been  made  on  the  projects  reported  in  the  letter  of  May 
12.  This  work  is  described  in  seme  detail  in  the  enclosed  expository  paper 
"Coupled  Oscillators  and  the  Design  of  Central  Pattern  Generators"  [1].  There 
are  a  large  number  of  related  projects  that  are  being  worked  on 
simultaneously;  the  main  reason  for  this  is  that  interactions  with  the 
experimentalists  have  raised  questions  relevant  to  the  design  and 
interpretation  of  experiments,  and  it  was  important  to  investigate  enough  of 
each  to  allow  the  experiments  to  go  on.  The  projects  and  their  current  state 
of  completion  is  as  follows: 


Completed  :  Two  expository  papers  (enclosed),  written  for  different  audiences. 

Structure  and  Function  in  an  Oscillating  Neural  Network"  is  to  appear  in 
Computational  Neurosciences.  E.  Schwartz,  ed.  MIT  Press,  1988.  Reference  [1] 
is  to  appear  in  Math  Biosciences  and  will  be  reprinted  in  Nonlinearity  in 
Biology  and  Medicine.  North  Holland  Publishing  Company,  1988 


Close  to  Completion:  (available  in  draft,  but  not  final,  form) 

1.  "Phase  transitions  and  other  phenomena  in  chains  of  oscillators". 

This  paper  explores  generalizations  of  the  work  in  our  paper  "Symmetry  and 
phaselocking  in  chains  of  weakly  coupled  oscillators"  [2] ,  on  which  most  of 
our  work  on  central  pattern  generators  is  based.  The  generalizations  were 
motivated  by  the  design  of  feasible  experiments,  many  of  which  have  recently 
been  done,  or  are  scheduled  to  be  performed  in  the  near  future.  (See  Section 
VI.)  The  generalizations  deal  with  the  effects  of  local  regions  of  higher 
frequency,  and  coupling  that  is  not  the  same  everywhere  along  the  cord.  It 
also  explores  the  implications  of  coupling  for  which  the  downward  (caudal) 
coupling  is  quite  different  from  the  upward  (rostral)  coupling.  New 
mathematical  phenomena  have  been  identified,  notably  the  existence  of 
solutions  that  exhibit  a  phase  transition  as  some  parameter  is  varied.  (Here 
a  phase  transition  refers  to  a  large  change  in  the  qualitative  behavior  of  the 
solution  that  occurs  over  a  tiny  change  in  the  parameter  if  the  number  of 
oscillators  is  large. )  These  results  yield  predictions  about  the  split  bath 
experiments  (Section  VI)  as  well  as  giving  a  framework  within  which  detailed 
results  of  those  experiments  can  yield  information  about  differences  between 
rostral  and  caudal  coupling,  and  differences  in  natural  frequencies  along  the 
cord. 


2.  "Averaged  and  unaveraged  systems  of  coupled  oscillators" 

This  paper  has  two  major  themes.  The  first  is  that  the  form  of  interaction  of 
biological  oscillators  commonly  used  in  models  of  forced  oscillators,  when 
applied  to  reciprocally  coupled  systems,  leads  to  .destruction  of  the 
rhythmicity.  This  is  very  robust,  and  holds  for  large  classes  of  models, 
including  standard  models  [3,4].  Particularly  susceptible  to  such  "oscillator 
death"  are  oscillators  interacting  by  means  of  a  single  large  pulse  in  each 
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oscillator  cycle.  The  second  these  is  that  this  undesirable  outcome  can  be 
avoided  if  the  interactions  take  place  in  a  sore  gradual  Banner,  spread  out 
within  the  cycle,  and  we  discuss  neural  networks  that  nay  be  doing  this. 
Technically,  we  show  that  if  the  interactions  cure  spread  out,  as  above,  then 
the  systea  behaves  like  one  in  which  the  coupling  is  averaged  over  the  cycle; 
we  show  that  for  such  systeas,  "oscillator  death"  cannot  occur.  This  work 
serves  to  support  the  assuaptions  aade  in  [2],  in  which  the  aethod  of 
averaging  was  used  to  reduce  a  seeaingly  unattackable  aatheaatics  problem  to 
one  for  which  analytic  conclusions  could  be  drawn. 

3.  With  D.  Aronson,  "Aaplitude  response  of  coupled  oscillators" 

This  paper  is  an  investigation  of  coupled  oscillators,  each  near  a  Hopf 
bifurcation,  so  that  the  strength  of  coupling  is  comparable  to  the  strength  of 
the  attraction  to  the  liait  cycle.  This  paper  has  been  substantially 
increased  and  iaproved  since  the  last  report.  Among  the  results  are  the 
existence  of  paraaeter  regimes  for  which  the  origin  is  stabilized  by  the 
interaction,  and  is  the  only  stable  solution.  We  are  thus  able  to  give  an 
analytical  explanation  for  a  phenomena  that  had  been  found  numerically  in 
several  models  of  cheaicals  oscillators.  We  are  also  able  to  deaonstrate 
bistability  of  a  phaselocked  solution  and  an  unlocked  "drift"  solution,  and 
bistability  of  a  phaselocked  solution  and  a  critical  point,  phenomena  which 
have  not  previously  been  found  analytically  in  coupled  oscillators. 


In  progress: 

1.  With  W.  Zhang  (graduate  student),  "Multiple  coupling  in  chains  of 
oscillators".  Neurons  froa  a  given  segment  of  the  lamprey  spinal  cord  have 
processes  that  extend  for  many  segments.  We  have  modelled  this  by  considering 
chains  of  oscillators  in  which  each  oscillator  is  connected  to  more  than  one 
of  its  neighbors,  as  in  the  leech  central  pattern  generator  [5].  We  show  that 
the  conclusions  of  [2]  still  hold  for  multiple  coupling.  In  addition,  we  show 
that  the  multiple  coupling  serves  to  buffer  the  solutions  against  variation  in 
frequency  of  the  local  oscillators.  Finite  chains  with  multiple  coupling  also 
appear  to  exhibit  edge  effects.  These  theoretical  results  help  to  explain  the 
otherwise  puzzling  results  of  A.  Cohen  on  the  measured  frequencies  and  phase 
differences  in  the  isolated  spinal  lamprey  preparation.  (See  Section  V.l  and 
[1,  Section  40])  The  latter  theoretical  results  are  preliminary;  if  they  turn 
out  to  be  robust,  they  may  give  unexpectedly  explicit  predictions  about  the 
anatomy  of  the  CPG  network.  In  particular,  the  mathematics  now  suggests  that 
the  position  of  the  "dips”  in  phase  differences  found  by  Cohen  are  a  measure 
of  how  aany  segments  are  connected  to  each  local  oscillator.  More  theoretical 
and  experimental  work  is  planned. 

2.  Work  on  "Periodically  forced  chains  of  oscillators"  was  motivated  by  a 
recent  paper  of  K.  Sigvardt  and  A.  McClellan  [6],  who  are  now  both  members  of 
the  Center.  Sigvardt  and  McClellan  investigated  sensory  feedback  in  the 
isolated  lamprey  spinal  cord  preparation,  using  vhe  aechanoreceptors  thought 
to  be  so-called  "edge  cells”  on  the  margins  of  the  spinal  cord.  They  were 
able  to  mechanically  force  a  few  segments  of  either  end  of  a  portion  of  a 
cord,  and  show  that  the  remaining  stationary  segments  could  be  entrained  by 


the  intersegmental  coordinating  system.  To  understand  the  details  of  their 
results,  we  modelled  the  situation  as  a  forced  chain  of  oscillators,  and  have 
preliminary  analytical  and  numerical  results.  We  found  these  results 
surprising;  they  assert  that  it  would  be  very  difficult  to  attain 
phaselocking,  as  was  seen  in  the  experiments,  for  a  range  of  forcing 
frequencies  and  with  forcing  at  either  the  rostral  or  the  caudal  end,  if  the 
rostral  and  caudal  coupling  are  qualitatively  the  same.  (See  Section  5C  of 
[1].)  The  analysis  shows  that  detailed  information  about  the  phase  lags  from 
forcing  experiments  have  the  potential  of  yielding  a  lot  of  information  about 
the  structure  of  the  C.P.G.  Because  of  this,  Sigvardt  and  T.  Williams 
(another  new  member  of  the  center)  have  undertaken  more  detailed  entrainment 
experiments.  (See  Section  VI.) 

3.  Work  on  "Nondiscrete  chains  of  oscillators"  is  motivated  by  the  need 
to  deal  with  the  fact  that  the  lamprey  and  other  vertebrate  spinal  cords  are 
not  composed  of  discrete  oscillating  units,  as  is  the  leech.  Up  until  now, 
that  has  been  outside  the  reach  of  any  modelling  efforts.  Our  work  is  based 
on  our  conjecture,  which  has  both  a  theoretical  and  experimental  basis,  that 
the  oscillating  units  are  composite.  (See  Section  5D  of  [2].)  Preliminary 
analytical  and  numerical  results  suggest  that  if  the  coupling  among  the 
neurons  within  an  oscillator  are  comparable  in  strength  to  coupling  among 
different  segments,  the  resulting  circuit  can  behave  like  a  chain  of  discrete 
oscillators,  while  not  anatomically  having  such  a  modular  structure. 

4.  A  collaboration  with  B.  Harder  has  begun.  We  are  interested  in  the 
differential  effects  of  standard  conductance  changes  and  voltage  dependent 
conductances  on  the  behavior  of  burster  neurons  and  networks  that  contain  such 
neurons.  (See  Section  IV.)  Preliminary  numerical  work  on  this  has  been 
carried  out  by  S.  Strogatz,  a  post-doctoral  fellow  associated  with  the  center. 
It  has  been  seen  that  existing  models  of  bursting  cells  do  not  respond 
appropriately  to  simulated  modulation  (e.g. ,  changes  in  baseline  and/or 
periodicity  are  incorrect);  other  models  are  being  explored.  The  numerical 
work  has  recently  been  taken  over  by  I.  Epstein  (Department  of  Chemistry, 
Brandeis  University),  who  is  spending  his  year  as  a  Guggenheim  Fellow  visiting 
the  Center. 

5.  Work  has  begun  on  mathematical  aspects  of  the  interaction  between  the 
neuronal  control  and  the  mechanics  of  swimming  in  lamprey.  This  work  involves 
a  collaboration  with  T.  Williams,  a  new  member  of  the  Center,  who  is 
knowledgeable  both  about  the  lamprey  [7]  and  about  physics.  (See  Section  VI.) 
The  mathematical  work  is  now  being  done  mainly  by  S.  Laederich,  a  graduate 
student  at  Boston  University,  under  the  supervision  of  N.  Kopell. 

6.  S.  Strogatz,  a  postoctoral  fellow  with  N.  Kopell,  has  been  working  on 
questions  involving  oscillators  with  stochasticity.  With  R.  Mirollo,  he  has 
written  a  paper,  "Phase-locking  of  oscillators  with  random  intrinsic 
frequencies",  submitted  to  Physica  D. 
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IV.  Work  from  the  Marder  Lab 

OVERVIEW 

We  are  interested  in  understanding  how  rhythmic  motor 
patterns  are  produced  by  the  nervous  system.  We  want  to 
define  the  mechanisms  by  which  rhythmic  motor  patterns  are 
modulated  by  neural  inputs,  sensory  information  and 
hormones.  We  want  to  know  how  the  activity  of  single  cell 
neural  oscillators  is  modulated  by  neurotransmitters  and 
modulators.  To  address  these  questions  we  use  the  pyloric 
rhythm  of  the  crustacean  stomatogastric  ganglion  as  a  model 
preparation. 

1.  IN  VIVO  RECORDINGS  OF  ONGOING  PYLORIC  ACTIVITY 

One  of  the  aims  of  the  initial  proposal  was  to  obtain 
good  quality  recordings  of  the  pyloric  rhythm  in  alive  and 
behaving  animals.  These  recordings  are  intended  to  provide 
important  baseline  data  about  the  frequency  and  variability 
of  the  pyloric  rhythm  under  a  variety  of  physiological 
conditions.  We  have  successfully  started  these  experiments, 
and  happily  find  that  the  crab  preparation  is  eaily  amenable 
to  the  implantation  of  chronic  recording  electrodes,  and 
that  data  are  reasonably  easy  to  obtain  in  these  in  vivo 
recording  configurations. 

Recordings  are  made  with  teflon  coated  silver  wires 
inserted  through  the  carapace  at  known  sites.  (Control 
experiments  were  performed  to  determine  where  the  muscle 
insertions  are  located.)  Visual  control  is  used  for 
electrode  insertion,  and  recordings  are  made  during  the 
insertion  process  to  enable  us  to  judge  when  the  electrode 
has  made  contact  with  the  appropriate  muscle.  When  the 
electrode  is  properly  positioned,  the  electrode  is  fixed  in 
place  with  surgical  cement,  and  the  animal  can  be  allowed  to 
move  about.  Injections  into  the  animal  are  made  directly 
into  the  circulation  in  the  region  of  the  heart. 

Figure  1  shows  an  example  of  recordings  made  from  two 
muscles  before  and  after  injection  of  proct.olin.  These 
recordings  clearly  indicate  that  proctolin  injection  into 
the  intact  and  behaving  animal  produces  dramatic  increases 
in  the  frequency  of  the  pyloric  rhythm,  as  well  as  also 
turning  on  activity  in  the  gastric  rhythm.  Figure  2  shows  a 
plot  of  the  effects  of  proctolin  injection  of  pyloric 
frequency  taken  from  records  such  as  those  in  Figure  1. 

These  experiments  are  in  progress  at  this  time.  These 
experiments  are  being  carried  out  by  Mr.  Robert  Zarum,  an 
undergraduate  who  is  currently  working  on  this  project  for 
his  Senior  Honors  Thesis  with  the  help  of  Dr.  Pierre 
Meyrand,  postdoctoral  fellow.  The  methodology  for  the 
collection  and  analysis  of  these  data  is  now  in  place. 


Figure  1. 

EMG  recordings  from  behaving  crab.  cpvl  muscle  recording 
shows  the  activity  of  the  pyloric  rhythm,  gml  muscle 
recording  shows  the  activity  of  the  gastric  mill  rhythm. 
Proctolin  injection  into  the  heart  at  t=0. 


Figure  2 . 


During  this  next  year  we  will  be  collecting  data  from  a 
large  number  of  animals,  analyzing  these  data,  and  preparing 
them  for  publication. 

2.  COUPLING  THE  GASTRIC  AND  PYLORIC  OSCILLATORS 

The  stomatogastric  ganglion  produces  two  different 
rhythmic  motor  patterns,  with  very  different  intrinsic 
frequencies.  The  pyloric  rhythm,  which  we  have  most 
intensively  studied,  has  a  period  of  about  one  second,  while 
the  gastric  rhythm  usually  displays  a  period  of  about  10 
seconds.  It  has  been  long  known  that  the  gastric  and  pyloric 
rhythms  interact  with  each  other  (Mulloney,  1977) .  This  is 
an  example  of  interactions  between  two  oscillators  of  very 
different  intrinsic  periods,  and  therefore  understanding  how 
these  oscillators  interact  was  of  interest. 

Recently  we  have  been  studying  the  interactions  between 
the  gastric  and  pyloric  rhythms  in  the  crab,  Cancer 
borealis .  In  this  species  we  find  that  a  number  of  the 
neurons  of  the  ganglion  appear  able  to  function  as  part  of 
both  gastric  and  pyloric  circuits.  Actually,  it  appears 
that  one  group  of  neurons  that  form  a  neuronal  "network"  can 
be  configured  into  different  neuronal  circuits  with  entirely 
different  functions.  These  data  are  derived  from  the  work 
of  Dr.  P.  Meyrand  and  Mr.  J.  Weimann,  and  are  currently 
being  prepared  for  publication. 

3.  MODULATION  OF  THE  AMPLITUDE  AND  FREQUENCY  OF  AN 
OSCILLATORY  NEURON 

Many  rhythmic  processes  in  the  nervous  system  are 
thought  to  be  controlled  by  oscillatory  "bursting"  neurons 
(Harris-Warrick  and  Flamm,  1986) .  These  are  neurons  that 
rhythmically  depolarize,  produce  bursts  of  action  potentials 
and  then  repolarize.  A  great  deal  of  work  has  been  expended 
in  trying  to  understand  the  ionic  mechanisms  underlying 
these  endogenously  active  oscillatory  neurons.  More 
recently,  research  has  started  in  trying  to  understand  how 
the  activity  of  bursting  pacemaker  neurons  is  modulated  by 
neurotransmitters  and  hormones. 

The  pyloric  rhythm  of  the  stomatogastric  ganglion  has  a 
timing  pacemaker  neuron,  the  Anterior  Burster  (AB)  that  can 
produce  bursts  of  action  potentials,  even  when  isolated  from 
other  synaptic  inputs.  We  have  previously  shown  that  the 
isolated  AB  is  modulated  by  amines  and  peptides  (Marder, 
1987)  .  We  are  now  interested  in  comparing  the  action  of 
different  classes  of  modulatory  substances  on  the  AB  neuron. 
Specifically  we  have  become  interested  in  comparing  the 
action  of  a  substance  that  increases  the  conductance  of  the 
AB  neuron  as  would  a  classical  excitatory  transmitter  with  a 
substance  that  modulates  the  action  of  voltage-dependent 
conductances  or  that  acts  in  a  voltage-dependent  fashion. 
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Therefore  we  have  chosen  to  study  the  action  of  two 
cholinergic  agonists  on  isolated  AB  neurons.  The  first 
substance,  nicotine,  increases  the  conductance  of  the 
postsynaptic  membrane.  The  second  substance,  pilocarpine, 
is  a  muscarinic  agonist,  and  has  unconventional  properties, 
similar  to  those  of  muscarinic  agonists  in  many  vertebrate 
preparations  (Marder  and  Paupardin-Tritsch ,  1978) .  We  are 
interested  in  studying  the  effects  of  these  two  agents  on  a) 
burst  amplitude  b)  baseline  membrane  potential  c)  burst 
frequency  and  d)  burst  duration.  Figure  3  illustrates  the 
basic  difference  between  the  action  of  these  two  substances 
on  an  AB  neuron:  nicotine  depolarizes  and  excites  the  AB 
neuron  but  does  not  cause  it  to  produce  bursts,  while 
pilocarpine  elicits  bursting  behavior. 

Some  of  these  data  are  included  in  a  paper  currently 
being  written  (Marder  and  Meyrand,  1988) .  As  physiologists, 
we  are  interested  in  fully  understanding  how  bursting 
neurons  are  modulated,  and  therefore  these  data  are 
interesting  to  us  in  and  of  themselves.  However,  these  data 
are  also  extremely  important  to  another  part  of  this 
project,  a  collaboration  with  the  Kopell  group,  on  the 
mathematics  of  the  modulation  of  rhythmically  active 
bursting  neurons.  This  collaboration  has  already 
demonstrated  that  a  simple  mathematical  model  that  was 
thought  to  describe  adequately  the  behavior  of  an  bursting 
pacemaker  neuron  (Rinzel)  does  not  account  in  a  satisfactory 
way  for  the  action  of  modulatory  substances  (see  Section 
III)  . 

4.  AN  EXPERIMENTAL  AND  THEORETICAL  EXPLORATION  OF  THE 
DIFFERENCES  BETWEEN  DIFFERENT  KINDS  OF  MODULATORS  ON  NEURAL 
CIRCUITS. 

Work  on  the  lobster  and  crab  stomatogastric  nervous 
systems  has  shown  us  that  the  pyloric  rhythm  is  modulated  by 
at  least  ten  different  substances,  including  a  number  of 
amines  as  well  as  many  peptides  (Marder,  1987)  . 

Interestingly  enough,  many  of  these  substances  appear  to  act 
either  in  a  voltage-dependent  manner,  or  by  modulating 
voltage-dependent  conductances.  This  suggested  to  us  that 
there  may  be  a  fundamental  reason  why  modulation  of  a 
rhythmically  active  neural  circuit,  or  for  that  manner,  any 
neural  circuit  might  be  best  achieved  ,in  this  manner.  In 
fact,  we  have  found  that  when  substances  that  act  simply  to 
increase  the  conductance  of  one  or  more  of  the  neurons  the 
pyloric  network  are  applied  to  the  preparation,  then  it  is 
possible  to  completely  disrupt  the  normal  pyloric  rhythm,  as 
illustrated  in  Figure  4.  Therefore,  our  working  hypothesis 
is  that  some  of  the  "modulators"  that  act  at  the  level  of 
the  voltge  dependent  conductances  have  conditional  actions 
just  so  that  they  do  not  risk  disrupting  normal  rhythmic 
activity . 
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Figure  3. 


The  effects  of  bath  application  of  10~4  M  nicotine  and 
pilocarpine  on  the  isolated  AB  neuron.  Compare  the 
depolarization  produced  by  the  nicotine  with  the  induction 
of  oscillations  produced  by  pilocarpine. 
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It  is  our  hope  to  develop  these  ideas  more  fully.  We 
will  then  test  them  experimentally  as  well  as  is  possible. 
When  the  model  of  the  bursting  neuron  responds  correctly  to 
the  application  of  modulatory  substances  we  plan  to 
collaborate  with  the  Kopell  group  to  try  to  develop  a  formal 
understanding  of  the  limits  of  modulation  of  a  rhythmically 
active  network  by  different  classes  of  modulatory 
substances . 
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V.  Work  accomplished  in  the  Cohen  Laboratory 

1 .  Structure  and  organization  of  the  central  pattern 
generator  for  locomotion 

A  robust  conclusion  of  the  mathematical  modeling 
studies  has  been  the  fact  that  it  is  difficult  to  obtain 
uniform  phase  lags  along  a  chain  of  oscillators  whose 
natural  frequencies  differ  to  any  substantial  degree.  Cohen 
had  shown  that  the  frequencies  of  the  segmental  oscillators 
of  the  lamprey  were  typically  different  and  could  be  very 
different.  The  precise  nature  of  the  differences  was 
unclear,  but  the  presence  of  the  differences  was  not.  On 
the  other  hand,  it  has  been  part  of  the  "lore"  of  the 
isolated  lamprey  spinal  cord  preparation  that  the  phase 
delays  along  the  cord  are  constant  over  a  wide  range  of 
frequencies  of  bursting  as  well  as  uniform  along  the  cord. 
However,  the  evidence  for  the  latter  was  inadequate  either 
to  confirm  or  deny  the  contention.  Experiments  were 
undertaken  by  Cohen  to  examine  both  the  nature  of  the 
frequencies  of  the  segmental  oscillators  along  the  cord  and 
the  precise  pattern  of  phase  delays  among  pairs  of 
oscillators  at  each  level  of  the  cord. 

A.  Methods 

a.  Phase  lag  determinations 

Fifty  segment  cord  pieces  were  dissected  out  and 
placed  in  the  bath  with  6  recording  electrodes  to  monitor 
ventral  root  activity.  Recordings  were  made  relative  to  2 
reference  electrodes,  one  on  each  side  of  the  cord.  The  4 
test  electrodes  were  moved  along  the  cord  in  steps  of  3 
segments  and  25  cycles  recorded  in  each  position.  The  phase 
lags  were  measured  for  each  position  relative  to  the 
reference  roots  (fig.l).  The  successive  phase  lags  were 
subtracted  one  from  the  next  and  the  difference  divided  by 
the  number  of  segments  intervening.  This  was  taken  to  be 
the  phase  lag  per  segment  over  the  segments  between  the  2 
recording  sites.  The  computed  phase  lag  per  segment  was 
plotted  against  the  segment  number. 

b.  Oscillator  frequencies 

In  the  same  cords,  the  preferred  frequencies  of  the 
regional  groups  of  oscillators  was  determined  by  two 
different  methods,  1)  selective  activation  by  means  of  cord 
partition  with  vaseline  seals  and  2)  cutting  the  cord  in  2 
and  then  4  pieces. 

B.  Results 

<1.  Phase  lags  along  the  cord 

By  the  method  described  in  Sectior  A  it  was  found  that, 
in  fact,  the  phase  delays  along  the  cord  were  not  uniform 
(fig.  2A) .  The  figure  usually  given  for  the  phase  delays 
has  been  an  average  of  1%  of  the  cycle  per  segment  with  no 
systematic  variation  in  any  particular  region  of  the  cord. 


The  delays  actually  found,  although  quite  variable,  had  the 
following  pattern.  The  delays  at  both  edges  of  the  pieces 
were  near  1%,  but  often  rose  a  few  segments  farther  from  the 
edges.  Near  the  segments  10  from  the  edges  there  was  a 
large  dip  in  the  delays.  Indeed,  the  phase  lags  in  this 
region  were  often  negative.  Still  farther  in,  the  delays 
rose  above  1%  briefly  and  then  fell  back  to  the  expected 
value. 

To  rule  out  the  possibility  that  this  pattern  is 
related  to  anatomical  anomolies  in  the  50  segment  cord 
region  tested,  Cohen  then  did  the  same  type  of  experiment  on 
cord  pieces  dissected  out  10  segments  farther  caudally  than 
in  the  previous  experiments.  Although  dislocated  relative 
to  the  previous  cords,  the  second  set  of  cord  pieces  gave 
essentially  the  same  pattern  of  dips  and  rises  in  roughly 
the  same  segment  numbers  (fig.  2B) .  Thus,  there  appears  to 
be  some  "edge  effects"  in  those  segments  around  10  from  the 
end  of  a  given  piece  of  cord. 


h.  Oscillator  f recruencies 


Cohen  had  suggested  on  the  basis  of  earlier  preliminary 
results  that  there  was  a  monotonic  gradient  of  frequencies 
along  the  cord.  However,  she  has  now  found  that  there  can 


be  monotonic  gradients,  but  they  are  the  exception  rather 
than  the  rule;  more  common  is  for  the  periods  to  be  longer 
in  the  end  regions  and  steadily  rising  or  falling  in  the 
mid-regions  (fig.  3)  For  reasons  that  are  still  unclear  the 
two  methods  of  frequency  determination  don't  always  produce 
the  same  pattern  of  results.  Thus,  although  there  are 
differences  in  the  frequencies  of  the  oscillators  along  the 
cord,  the  magnitude  and  pattern  of  the  differences  will  only 
be  known  after  more  work,  including  work  to  clarify  the 
method  best  for  making  the  determination. 

This  work  has  been  reported  in  an  abstract  of  work  to 
be  presented  at  the  American  Neurosciences  meeting  in 
November,  1987  and  a  meeting  in  Florida  (NAME) ,  October, 
1987.  The  latter  will  be  accompanied  by  a  paper  to  be 
published  as  part  of  the  preceedings  of  the  meeting,  and 
will  include  an  Appendix  by  Kopell. 


C.  Future  work  on  these  questions 

In  new  modeling  work  precipitated  by  the  phase  lag 
data,  a  pattern  similar  to  the  one  described  has  been  found 
to  occur  with  multiple  neighbor  coupling.  In  the 
simulations  done  to  date  on  the  new  models,  in  some 
parameter  ranges  the  number  of  segments  at  which  the 
negative  dips  occur  corresponds  to  the  number  of  segments 
connected  on  each  side  to  a  segmental  oscillator.  Thus, 
these  models  predict  tentatively  that  there  are  coordinating 
fibers  connecting  up  to  10  contiguous  segments.  Becuase  of 
the  clearness  and  usefulness  of  this  prediction,  it  has  been 
decided  to  continue  this  line  of  investigation  in  greater 
detail . 


Nicholas  Mellen,  Cohen's  graduate  student,  will  do  this 
as  a  portion  of  his  graduate  research.  He  has  designed  and 
is  currently  building  a  set-up  which  will  include  electrode 
holders  for  easy  positioning  on  nearby  segments.  This  will 
allow  measurement  of  the  phase  lag  over  a  few  segments 
directly  instead  of  relative  to  a  reference  root  at  a 
distance;  and  may  thereby  reduce  the  variability  in  the  data 
as  well  as  speed  up  the  data  collection  nd  analysis.  The 
set-up  will  also  be  built  with  improved  partitions  for  the 
frequency  experiments.  It  is  hoped  that  this  will  improve 
these  data  and  facilitate  the  experiments. 

First  he  will  test  whether  the  effect  still  exists  on 
shorter  cord  pieces  (40,  30,  20  and  10  segments  long) .  For 
long  enough  pieces  the  mathematics  suggests  that  the 
negative  dips  should  remain  10  segments  from  the  end; 
numeric  simulations  will  be  done  on  short  pieces  of  cord  to 
see  if  the  model  matches  resultant  data.  If  the  modeling 
holds  up  in  this  regard,  then  this  type  of  experiment  will 
be  used  to  define  the  regional  anatomy  of  the  coordinating 
system.  For  example,  the  experiments  described  above  were 
all  done  on  rostral  pieces  of  cord.  It  is  unknown  whether 
the  more  caudal  portions  of  the  cord  have  the  negative  dips. 
If  they  don't  then  it  will  suggst  that  they  do  not  have  the 
same  anatomy  as  the  more  rostral  portions  used  in  the 
experiments  described  above. 

2 .  Regeneration 

A.  Larval  lampreys 

Studies  on  larval  lampreys  were  undertaken  by  Cohen, 
Mackler  and  Selzer  (1985)  to  confirm  that  the  functional 
regeneration  was  independent  of  mechanical  entrainment. 

This  possibility  was  raised  by  Grillner  and  his  associates 
in  an  abstract  in  which  they  reported  that  the  neuromuscular 
blocker,  D-tubocurarine,  disrupted  the  coordination  observed 
during  NMDA  induced  fictive  swimming  by  formerly  transected 
larval  cords.  To  test  whether  or  not  they  were  correct, 
Cohen  repeated  her  original  experiments  with  the  addition  of 
curare  to  the  bathing  solution.  She  found  that  curare  had 
no  discernible  effect  on  the  coordination  maintained  by 
regenerated  fibers  (fig.  4) .  The  only  problems  that  could 
have  accounted  for  Grillner 's  negative  results  were  that  his 
group  used  higher  doses  of  NMDA  than  she  found  optimal  and 
that  they  probably  did  not  wait  long  enough  for  the 
coordination  to  stabilize  (it  can  take  as  long  as  45 
minutes).  Cohen's  curare  work  has  been  submitted  for 
publication  to  the  Journal  of  Neurobiology . 

B.  Adult  lamprevs 

There  has  also  been  success  in  keeping  the  spinal  cord 
lesioned  adult  lampreys  alive.  They  have  now  been  surviving 
in  adequate  numbers  (>25%)  for  10  months.  The  critical 
factor  appears  to  be  feeding  the  animals  and  employing 
partial  lesions  instead  of  complete  transections .  The 
partial  lesions  allow  the  animals  to  feed  successfully  since 


they  are  predatory  on  free-swimming  fish  and  the  partial 
lesions  don't  interrupt  swimming  during  the  recovery  period. 
Most  laboratories  that  have  attempted  to  keep  the  animals  do 
not  have  the  facilities  to  feed  their  animals  and  the  young 
adult  lampreys  used  for  these  experiments  simply  cannot 
survive  without  food  on  a  regular  basis.  Cohen  expects  even 
greater  survival  in  the  coming  year  when  arrangements  will 
be  made  to  keep  the  prey  fish  available  throughout  the 
entire  required  recovery  period.  This  past  year  there  was  a 
hiatus  of  two  months  which  precipitated  the  deaths  of 
several  animals. 

Tests  of  the  lesioned  animals  have  now  begun  and  have 
yielded  interesting  results.  The  lesions  made  were  of  1) 
hemicords,  2)  lateral  tracts,  and  3)  medial  tracts.  In 
isolated  control  cords,  none  of  these  lesions  when  made 
acutely  disrupted  coordination  (cohen,  1987)  .  Regenerated 
fibers  of  only  the  medial  tracts  have  consistently  produced 
animals  (3  of  3)  with  firtive  swimming  coordinated  across 
the  lesioned  segment  after  the  intact  regions  were  cut 
acutely  in  the  bath.  In  all  cases,  the  bursting  was  never 
well  coordinated.  Typically,  it  would  lock  for  sometime, 
then  drift  apart  again  (fig.  5).  Based  on  modeling  studies 
(Cohen  et  al.,  1982)  these  observations  are  consistent  with 
weak  coupling,  much  weaker  than  that  observed  in  larval 
cords. 

These  observations  are  in  contrast  to  preliminary 
experiments  on  hemisected  and  laterally  lesioned  cords  which 
lead  to  much  more  complex  results.  To  date,  with  2  animals 
tested,  the  lateral  lesions  produced  one  cord  with  no  clear 
evidence  of  coupling  (fig.  6),  but  which  did  produce  fictive 
swimming  on  each  side  of  the  lesion.  The  hemisected  and  the 
other  laterally  lesioned  cord  could  not  be  induced  to 
produce  any  patterned  bursting  in  the  segments  surrounding 
the  lesion  site  although  the  segments  at  some  distance 
rostral  and  caudal  to  the  region  of  the  lesion,  when  cut 
apart  from  the  area  adjacent  to  the  lesioned  segments 
produced  stable  bursting.  The  fact  that  segments  at  a 
distance  would  produce  stable  bursting  while  those  segments 
in  the  vicinity  of  the  lesion  would  not  suggests  that  the 
regenerated  fibers  somehow  perturbed  the  central  pattern 
generator. 

C.  Future  work  on  regeneration 

These  observations  made  on  Petromvzon  marinus  are 
being  confirmed  with  the  remaining  animals.  They  are  also 
being  extended  to  include  animals  of  the  species 
Ichthvomvzon  unicuspis.  Preliminary  experiments  with  this 
latter  species  point  to  its  having  less  trouble  producing 
fictive  locomotion  regardless  of  the  lesion,  but  this  will 
be  tested  directly  in  2-3  months  when  they  have  had  adequate 
time  to  recover. 

The  observations  outlined  above  suggest  that  the  adult 
lamprey  spinal  cord  can  be  used  as  a  model  system  to  study 
the  factors  underlying  both  regeneration  that  restores 
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function  and  that  which  leads  to  disrupted  function. 

Possible  disruption  by  regeneration  has  been  played  down  by 
other  investigators  who  have  pointed  instead  to  the  apparent 
normalcy  of  the  regenerated  axons  and  their  connections 
(Mackler  and  Selzer) .  However,  the  comparison  between  the 
two  types  of  regeneration  will  be  exploited  by  Cohen's  lab. 
Tamara  Dobrof,  a  graduate  in  the  lab,  will  be  following  up 
on  these  experiments  both  anatomically  and  physiologically. 
She  will  learn  histological  techniques  (HRP  backfilling) 
necessary  for  examining  the  origin  and  destination  of  the 
regenerated  fibers.  Dr.  Karen  Sigvardt,  a  new  member  of  the 
Center,  has  agreed  to  teach  the  necessary  methods  to  Dobrof 
this  winter.  Since  it  is  intuitively  likely  that  the 
difference  between  the  lesion  results  lies  in  the  classes  of 
neurons  with  fibers  in  the  different  tracts  and  the 
connections  they  form  upon  regeneration,  the  initial  thrust 
of  her  work  will  be  to  classify  the  neurons  whose  fibers 
pass  in  the  different  fiber  tracts  in  control  and 
regenerated  cords. 

3.  Optical  recording  of  neuronal  activity  through  the  use 
of  voltage  sensitive  dyes 

The  equipment  necessary  for  the  optical  recording  of 
neuronal  activity  is  now  largely  assembled.  Unfortunately, 
the  analogue  to  digital  converter  being  built  by  an 
electronic  technician  at  Yale  University  is  still 
incomplete,  but  has  been  promised  in  3  months.  At  that 
time,  L.B.  Cohen,  Yale  University,  has  agreed  to  take 
responsibility  for  writing  the  programs  necessary  for  taking 
and  analyzing  data.  He  also  needs  the  equipment,  and  has 
agreed  to  share  all  developments  fully  with  us.  All  the 
other  equipment  has  been  successfully  brought  to  a  state  in 
which  recording  has  been  achieved  on  the  lamprey  spinal 
cord.  There  remain  improvements  to  be  made  in  the  recording 
arrangement,  but  it  is  now  possible  to  record  neuronal 
activity  optically  in  Cohen's  own  laboratory.  She  will 
begin  to  optimize  further  the  conditions  to  allow 
measurement  of  the  small  signals  produced  by  the  lamprey 
cells  post-synaptic  to  normal  and  regenerated  neurons.  Dr. 
Guan  Li  of  the  People’s  Republic  of  China,  will  join  her 
this  winter  to  help  in  this.  Dr.  Guan's  training  in 
neurophysiological  techniques  is  extensive  and  she  should  be 
able  to  contribute  to  the  work  with  a  minimum  of 
supervision. 

Dr.  Peter  Brodfuehrer  will  also  work  for  a  short  time 
with  Cohen  to  try  to  record  from  leech  neurons  on  the 
optical  set-up.  Brodfuehrer  did  his  thesis  on  leech  central 
pattern  generation  of  locomotion  and  will  try  to  extend  his 
studies  through  the  use  of  optical  recording  techniques.  He 
will  begin  by  screening  dyes  to  find  those  that  have  little 
or  no  effect  on  the  CPG.  The  initial  testing  will  be  done 
this  winter.  If  these  tests  are  successful,  he  and  Cohen 
will  submit  a  small  grant  application  to  support  his  work  in 
her  lab  on  the  leech  CPG. 
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FIGURE  LEGENDS : 


FIGURE  1 :  Diagram  of  method  for  obtaining  segment-by- segment 
phase  delays. 

FIGURE  2:  Segment-by-segment:  phase  delays  along  50  segment 
cord  pieces.  Upper  graph  is  for  segments  0-50  beginning  just 
behind  gills.  The  graphs  for  each  cord  vary  except  for  the 
two  rises  and  dips  centered  near  10  segments  from  each  end. 
The  arrow  indicated  the  mean  of  the  segment  numbers  that  each 
cord  has  a  negative  phase  delay  in  this  region.  The 
horizontal  bar  is  the  standard  deviation.  One  cord  had  no 
negative  dip  in  the  caudal  region.  Lower  graph  is  the  same 
format  but  for  cord  pieces  10  segments  caudal  to  those  used 
to  obtain  upper  graph.  Notice  the  arrows  are  still  about  10 
segments  from  the  ends,  even  though  the  cord  piece  is  now 
segments  10-60. 

FIGURE  3:  Cycle  Periods  for  cord  pieces.  The  mid-segment  of 
the  piece  is  plotted  against  cycle  period.  The  cords  are  the 
same  as  those  used  for  figure  2  data. 

FIGURE  4:  Bursting  by  isolated  larval  lamprey  cords  that  had 
been  transected.  The  upper  trace  is  a  segment  rostral  to  the 
lesion;  lower  trace  is  caudal  to  the  lesion.  Bursting  was 
induced  by  D-glutamate  plus  curare.  Lower  graph  is  a 
histogram  of  the  phases  which  clearly  center  at  a  stable 
value . 

FIGURE  5:  Phase  delays  over  time  of  an  isolated  adult  cord. 
The  bursting  was  by  segments  rostral  and  caudal  to  a  lesioned 
segment  in  which  only  regenerated  medical  tract  and  gray 
region  were  left  intact.  The  phase  lags  are  variable  but 
clearly  show  regulation  which  is  pulling  the  activity  back 
and  forth  to  some  preferred  value. 

FIGURE  6:  Same  as  figure  5  except  lesioned  segment  left  only 
lesioned  lateral  tracts  intact.  The  phase  delays  demonstrate 
almost  pure  drift  of  rostral  and  caudal  segments. 


PAPERS  SUBMITTED: 


Evidence  for  functional  regeneration  in  lamprey  spinal  cord  persists  in  the  absence 
of  reflex  movements.  To:  Journal  of  Neurobiology. 


ABSTRACTS  SUBMITTED: 


Myths  and  realities  of  frequency  gradients  and  constant  phase  lags  during  lamprey 
fictive  swimming.  Society  for  Neuroscience  Abstract. 

Evidence  for  functional  regeneration  by  adult  lamprey  spinal  cords.  VA-PVA  Symposium. 


PAPERS  IN  PREPARATOIN  OR  PLANNED: 


Studies  of  the  lamprey  central  pattern  generator  for  locomotion:  A  close  relation¬ 
ship  between  modeling  and  experimentation.  (Invited  for  meeting  on  Dynamic  Patterns 
in  Biology ) . 

Evidence  for  functional  regeneration  by  adult  lamprey  spinal  cords.  Science  (to 
be  submitted). 

Detailed  analysis  of  intersegmental  interactions  in  the  isolated  spinal  cord: 
Unexpected  systematic  variations  in  phase  relations.  Brain  Research  or  Neuroscience 
( to  be  submitted ) . 


V.  Work  accomplished  by  Sigvardt  and  Williams. 


1.  Control  and  generation  of  rhythmicity  in  the  spinal  cord .  (Alford, 
Sigvardt  and  Williams). 

In  previous  work  we  have  demonstrated  that  rhythmic  locomotory  output 
from  the  spinal  cord  is  elicited  by  N-methyl-D-aspartate  (NMDA)  receptor 
activation  (Grillner,  et  al.,  1981).  We  then  showed  that  there  are  neurons  in 
the  spinal  cord  that  have  NMDA- dependent  membrane  potential  oscillations 
(Sigvardt  et  al . ,  1985).  It  has  been  hypothesized  that  inhibitory  glycinergic 
transmission  plays  an  important  role  in  the  generation  of  rhythmicity 
providing  the  inhibitory  phase  of  the  membrane  potential  oscillations  in 
ventral  horn  neurons.  However,  this  year  Alford  and  Williams  have 
demonstrated  that  the  spinal  cord  networks  produce  a  robust  rhythmic  output 
when  glycinergic  transmission  is  blocked  by  strychnine  and  hypothesize  that 
burst  termination  is  the  result  of  voltage -dependent  processes  inherent  to 
neurons  within  the  network.  Glycinergic  transmission  would  then  provide  a 
timing  signal  to  ensure  that  the  networks  on  the  opposite  sides  of  the  spinal 
cord  are  activated  in  antiphase.  We  are  now  doing  experiments  to  verify  the 
mechanism  for  burst  termination  in  the  presence  of  strychnine.  With 
glycinergic  transmission  blocked,  the  activity  on  opposite  sides  of  the  spinal 
cord  is  synchronous.  In  the  coming  year,  we  have  planned  experiments  to 
determine  the  locus  of  this  coupling  and  the  neurons  that  might  provide  the 
glycinergic  timing  signals  to  the  oscillators  during  naturally-  induced 
activity. 


2.  Forcing  experiments :  forced  oscillators  and  intersegmental  coupling. 
(Sigvardt  and  Williams) . 

In  these  experiments  a  movement -signal  was  applied  to  one  end  of  a  piece 
of  spinal  cord/notochord  containing  a  defined  number  of  segments  from  a 
defined  anatomical  location  within  the  cord  and  the  range  of  forcing 
frequencies  over  which  the  spinal  cord  oscillators  can  be  entrained  was 
determined  and  the  intersegmental  phase  lags  measured.  It  was  critical  for 
these  experiments  that  the  applied  movements  were  carefully  controlled.  The 
movement  was  produced  by  a  Gould  pen  motor  supplied  by  the  center  and 
controlled  by  a  system  designed  by  Dr.  McClellan  and  me.  The  usefulness  of 
these  experiments  for  understanding  the  organization  of  the  spinal  locomotory 
networks  depends  on  the  previous  observation  of  McClellan  and  Sigvardt  that 
mechanoreceptors  in  the  spinal  cord  transmit  movement -related  feedback  signals 
onto  local  oscillators  and  the  changes  induced  at  the  local  level  are 
transmitted  through  the  intersegmental  coordinating  system.  Therefore,  these 
forcing  experiments  provide  information  about  both  the  strength  of  the 
feedback  system  and  the  intersegmental  coordination  system.  The  experiments 
were  designed  to  address  the  following  questions:  1).  What  factors  influence 
the  range  of  frequencies  over  which  fictive  locomotion  can  be  entrained?  a) 
number  of  spinal  cord  segments?  b)  rostral  versus  caudal  forcing?  c) 
anatomical  location  (rostral  segments  versus  caudal  segments)?  2).  Does 
forcing  affect  intersegmental  coordination  (phase  lags)?  a)  are  the  effects 
local  or  global?  b)  do  these  changes  in  intersegmental  coordination  differ 
with  rostral  as  compared  to  caudal  forcing?  3) .  What  are  the  phase 
relationships  between  movement  and  the  ventral  root  bursts?  a)  Are  these 
relationships  related  to  any  of  the  above  factors;  e.g.  does  the  relationship 


between  the  movement  phase  and  forcing  frequency  affect  the  entrainment  range? 
b)  How  do  these  relationships  compare  to  those  in  the  intact  animal  when  the 
feedback  loop  is  unbroken?  The  experiments  were  completed  this  summer  and  we 
are  now  analyzing  the  data.  Plans  for  further  experiments  of  this  type  depend 
on  the  results  of  this  analysis  which  will  be  interpreted  within  the  context 
of  the  coupled  oscillator  theories  of  Kopell  and  Ermentrout. 


3.  Split-bath  experiments  (Sigvardt,  Williams  and  Kopell). 

The  coupled  oscillator  theory  developed  by  Kopell  and  Ermentrout  allows 
one  to  make  predictions  about  the  effects  of  various  experimental 
manipulations  on  the  intersegmental  coupling  within  the  spinal  cord  networks. 
I  have  done  a  series  of  experiments  in  which  intersegmental  coupling  is 
examined  when  the  intrinsic  frequency  of  the  oscillators  in  one  half  of  the 
spinal  cord  piece  is  changed  relative  to  that  in  the  other  half  by  bathing  the 
two  halfs  of  the  cord  in  different  concentrations  of  excitatory  amino  acid. 
If  the  inherent  frequency  of  the  rostral  half  of  the  spinal  cord  is  higher 
than  the  caudal  half,  the  phase  lags  along  the  cord  are  longer  than  normal 
whereas  if  the  inherent  frequency  of  the  caudal  half  is  higher  than  the 
rostral  half,  there  is  a  negative  phase  lag  (backward  swimming)  at  the 
transition  between  caudal  and  rostral  and  phase  lags  more  rostrally  are 
reduced  to  zero.  We  are  also  planning  experiments  to  raise  the  inherent 
frequency  in  a  few  segments  in  the  middle  part  of  the  cord  to  see  if  the  waves 
will  propagate  outward  from  this  point  as  the  theory  predicts.  The  results  of 
these  experiments  should  provide  information  about  inherent  frequency  and  the 
strength  and  timing  of  ascending  and  descending  coupling  in  the  maintenance  of 
intersegmental  coordination. 


4.  The  relationship  between  muscle  activation  and  movement .  (Williams). 

In  locomotion  the  movement  produced  is  the  result  of  the  pattern  of 
muscle  activation  generated  by  the  central  nervous  system.  Williams  in 
collaboration  with  Sten  Grillner  and  his  colleagues  has  examined  the  match 
between  muscle  activation  and  movement  in  the  lamprey  and  the  trout,  two  fish 
that  have  very  different  forms  of  swimming.  They  have  found  that  the 
activation  of  the  contractile  elements  in  both  animals  is  approximately  in 
phase  with  the  velocity  of  that  part  of  the  body  of  the  fish  that  is 
performing  the  most  work  against  the  water.  Dr.  Williams  has  developed  a 
theory  relating  muscle  activation  to  movement  by  considering  the  body  of  the 
animal  as  a  damped  oscillator  with  its  own  natural  frequency  and  then  viewing 
the  tension  developed  by  the  muscle  activation  as  a  driving  force  extrinsic  to 
the  passive  oscillator.  In  this  way  the  phase  coupling  between  the  forcing 
function  (muscle  activation)  and  the  induced  oscillations  (movement)  at 
different  positions  along  the  body  can  be  compared  to  the  phase  coupling  in  a 
damped  linear  oscillator  being  driven  by  a  simple  periodic  force  at  different 
frequencies.  This  theory  provides  an  explanation  of  why  the  phase  between 
muscle  activation  and  movement  changes  along  the  body  of  the  animal. 
Mathematical  questions  raised  by  this  conjecture  are  besing  investigated  by  N. 
Kopell  and  S.  Laederich. 


Studies  of  muscle  physiology:  the  relationship  between  muscle  activation 
and  force  generation ,  (Curtin,  Sigvardt  and  Williams). 

In  order  to  fully  understand  the  effects  of  the  changing  phase  between 
muscle  activation  and  movement  down  the  length  of  the  body  one  need  to  know 
the  relationship  between  muscle  activation,  tension,  shortening  and  time. 
Williams  and  Sigvardt  will  work  collaboratively  this  year  with  Curtin  to 
describe  these  relationships  in  lamprey  muscle  so  that  future  models  of  the 
neural  basis  of  movement  can  include  these  parameters. 

6.  Methods  of  data  inspection,  display  and  analysis. 

a.  Strip  chart  recorder:  All  of  the  experiments  done  in  my  laboratory 
involve  recording  extracellularly  from  between  2  and  six  ventral  roots  along 
the  length  of  the  spinal  cord  in  the  in  vitro  lamprey  spinal  cord  preparation 
during  fictive  locomotion  to  test  the  effects  of  various  pharmacological 
manipulations,  anatomical  lesions  and  applied  movements  on  the  intersegmental 
coordination  along  the  length  of  the  spinal  cord.  I  had  in  my  laboratory  all 
the  equipment  necessary  to  do  the  experiments  except  a  chart  recorder  for  data 
display  and  analysis.  Since  the  value  of  the  these  experiments  depends  on 
precise  quantitative  analysis  of  intersegmental  phase  lags  over  many  cycles  at 
at  least  two  and  as  many  as  six  points  along  the  length  of  the  spinal  cord,  a 
hard  copy  of  all  of  the  data  was  needed  from  which  to  make  measurements  of 
phase  lags.  The  center  purchased  a  refurbished  Gould  260  chart  recorder  and 
accompanying  pen  motors  and  amplifiers  from  Pacer  Scientific  Instruments.  The 
hard  copy  of  the  data  can  be  made  on-line  or  off-line  from  my  FM  tape  recorder 
onto  this  recorder. 

b.  ASYST  training  course .  (Sigvardt). 

On  September  15-17,  1987  I  attended  a  training  seminar  to  learn  ASYST. 
ASYST  is  a  scientific  software  system  that  allows  acquisition,  analysis  and 
presentation  of  data.  Furthermore,  ASYST  is  a  programming  language  that  is 
extremely  easy  to  use  to  develop  simple  commands  or  complex  programs  for  data 
acquisition  and  analysis.  The  acquisition  system  allows  one  to  collect  data 
from  most  of  the  widely  used  A/D  boards  as  well  as  to  create  data  files  from 
other  systems  that  have  been  stored  on  floppy  discs.  The  analysis  package 
allows  one  to  edit  data  files  and  easily  do  any  type  of  mathematical  analysis 
from  simple  arithmetic  to  Fourier  analysis,  statistics  and  curve  fitting.  The 
graphics  package  provides  for  simple  plots  as  well  as  contour  plots  and 
several  types  of  charts. 

My  attendance  at  this  course  was  sponsored  by  the  center.  I  am  using 
ASYST  for  data  analysis  in  my  laboratory  and  plan  to  teach  the  system  to  the 
other  members  of  the  group.  In  all  of  our  experiments  we  collect  data  on 
magnetic  tape.  Then  relevant  segments  of  the  data  can  be  transferred  to 
floppy  disc  for  data  analysis  using  the  ASYST  system.  Since  so  much  of  the 
data  collected  is  of  interest  to  all  the  members  of  the  group,  we  think  that 
it  is  important  to  standardize  methods  of  data  analysis  and  display.  With 
ASYST  we  will  be  able  to  send  data  files  back  and  forth  between  laboratories 
on  floppy  disc  and  develop  programs  for  data  analysis  and  display  that  can  be 
used  by  all  of  us.  In  my  opinion  ASYST  is  an  excellent  system  because  it  has 
so  many  capabilities,  is  easy  to  use  and  the  staff  in  the  technical  support 
office  are  very  helpful. 
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VII-  Coordinating  activities  of  the  Center. 


1.  Conaolidation  and  Collaboration.  Through  the  efforta  of  the  Center,  aoat 
of  the  major  reaearchera  (west  of  Stockholm)  inveatigating  the  central  pattern 
generator  of  the  lamprey  are  working  together  in  a  joint  effort  to  understand 
thia  important  Model  system.  Three  auch  inveatigatora  have  joined  the  center 
in  the  paat  year:  K.  Sigvardt,  T.  Williams  and  A.  McClellan.  Many 
collaborationa  between  members  of  the  center  have  foraed  over  the  year,  with 
joint  work  performed  and  research  papers  planned.  These  include 

a.  Entrainment  experiaents.  Sigvardt  and  William  have  worked  together 
to  generate  new  data.  (See  Section  VI.)  A  joint  paper  with  Kopell  and 
Eraentrout  ia  planned  about  the  iaplications  of  the  data. 

b.  Split  bath  experiaents.  Sigvardt  has  per foraed  aoae  of  these 
experiaents,  and  ia  planning  others.  (See  Section  VI.)  A  joint  paper  with 
Kopell  and  Eraentrout  is  planned  about  the  theoretical  iaplications. 

c.  Kopell  plana  to  write  an  Appendix  to  a  paper  by  Cohen  discussing  the 
theoretical  iaplications  of  tue  edge  effects  found  in  her  data  on  phase  lags. 
(See  Sections  III  and  V.) 

d.  McClellan  has  agreed  to  provide  Williaaa  with  aoae  data  on  eechanical 
properties  of  intact,  spinal ized  and  dead  laapreys,  information  that  is 
iaportant  in  the  theories  of  Williaaa  and  her  mathematical  collaborators. 

e.  Marder  ia  working  with  Kopell  and  Epstein  on  questions  relating  to  the 
effects  of  neuroaodulators  on  the  output  of  the  lobster  stoaatogastric  central 
pattern  generator.  (See  Sections  III  and  IV.) 


2.  Training  of  graduate  students,  postdoctoral  associates  and  Center 
personnel  was  shared  by  the  aeabers  of  the  Center.  For  exaaple 

a.  S.  Alfords,  a  graduate  student  of  Williams,  is  spending  time  in  the 
laboratory  of  Sigvardt,  working  on  a  new  preparation  which  will  avoid  some  of 
the  technical  problems  of  working  with  the  isolated  spinal  cord.  So  is  M. 
Dobrof,  a  student  of  Cohen,  who  will  learn  some  techniques  in  histology. 

b.  The  Center  sponsored  the  attendance  of  Sigvard  at  a  software  training 
course,  in  which  she  learned  techniques  valuable  to  other  members  of  the 
Center,  and  which  she  expects  to  teach  to  other  members. 

c.  S.  Laederich,  a  graduate  student  of  Kopell,  is  working  on  mathematical 
issues  raised  by  the  theories  of  Williams. 


3.  The  funding  of  equipment,  training  and  travel  costs  of  new  members  of  the 
Center  was  accomplished  by  reallocation  of  some  of  the  money  in  the  core 
budget . 


4.  All  but  one  member  of  the  center  convened  from, Aug.  25  to  Aug.  28  to  take 
part  in  an  informal  workshop.  This  workshop  functioned  for  the  Center  as  a 
large  laboratory  meeting,  in  which  recent  data  and  unanswered  questions  were 
discussed  in  detail.  Some  of  the  themes  of  that  meeting  are  in  the  Informal 
Agenda,  in  the  Appendix. 


5.  The  Center  sponsored  a  widely  publicized  and  well-attended  sesinar  series 
on  "Dynaaical  Systems  in  Biology.”  Speakers  were  S.  Strogatz  (Boston 
University),  N.Kopell  (Boston  University),  J.  Rinzel  (N.I.H.),  D.  Tank  (A.T.&T 
Bell  Labs),  S.  Geaan  (Brown  University),  S.  Grossberg  (Boston  University), 

G.B.  Eraentrout  (University  of  Pittsburgh),  C.  Peskin  (Courant  Institute),  A. 
Perelson  (Los  Alaaos)  and  T.  Sejnowski  (Johns  Hopkins).  This  series 
eaphasized  oscillatory  phenoaena  and  network  phenomena  (neural  networks  and 
iaaunological  networks) . 


6.  The  Center  hosted  aany  visitors,  in  addition  to  the  ones  listed  above. 
Visitors  this  acadeaic  year  included  D.  Rue lie  (Institute  dec  Hautes  Etudes 
Scientifiques) ,  D.  Aronson  (University  of  Wisconson) ,  P.  Getting  (University 
of  Iowa),  T.  Williaas  (St.  Georges  Medical  School), K.  Sigvard,  (University  of 
California,  Davis),  A.  Stretton  (University  of  Wisconsin),  R.  Bisenberg  (Rush 
Medical  Center)  and  H.  Maturana  (University  of  Santiago).  C.  Rovainen 
(Washington  University  Medical  School)  spent  his  Fall  sabbatical  in  Boston, 
and  was  an  adjunct  aei&er  of  the  Center.  Another  such  adjunct  aeaber,  who 
partipates  in  sost  center  activities,  is  S.  Giszter,  (Bizzi  Lab,  M.I.T.).  Soae 
of  these  visits  were  the  occasion  of  the  gathering  of  outlying  aeabers  of  the 
Center  for  sini-workshops . 


Appendix 


Center  Workshop  -  Tues  Aug.  25  -  Fri  Aug.  28 
Informal  Agenda 


1.  Overview  of  lamprey  swinging  work.  What  can  the  reduced  preparation  tell 
us?  Problems  and  pitfalls  of  reduced  preparations.  Issues  of  central 
vs.  peripheral  control,  role  of  higher  structures.  Where  are  various 
kinds  of  regulation  performed?  Etc.  This  discussion  is  meant  to 
introduce  these  topics,  but  not  exhaust  them  (or  us!);  we  will  come  back 
to  them  as  necessary. 

2.  Some  relevant  mathematics.  An  informal  tutorial  on  some  mathematics  of 
oscillators,  concentrating  on  results  that  appear  to  be  directly  relevant 
to  experiments  that  have  been  done,  or  in  principle  can  be  done. 

3.  Detailed  discussion  of  experimental  results  on  the  isolated  spinal  cord. 
Interpretation  of  results.  Issues  involving  manipulation  of  data. 
Discussion  of  further  experiments  that  would  be  useful. 

A.  Structure  of  intersegmental  coordinating  system,  frequency  and  lesion 

experiments 

B.  Entrainment  experiments 

C.  Effects  of  changes  in  pharmacology  on  phase-lags 

D.  Intracellular  work,  discussion  of  effects  of  NMDA  and  strychnine 

4.  Initiation  of  motor  activity  at  different  levels. 

5.  Neuromodulation  of  small  neural  networks.  This  is  work  motivated  mostly 
by  work  on  invertebrates,  such  as  the  lobster  stomatogastric  ganglion. 
However,  principles  appear  to  be  emerging  that  are  likely  to  be  more 
generally  relevant. 

6.  Regeneration  experiments.  Central  vs.  peripheral  factors  in  behavioral 
recovery. 

7.  Biomechanics  and  muscle.  Neural  waves  and  mechanical  waves. 
Interpretation,  modelling,  discussion  of  further  experiments 


Participants:  N.  Kopell,  A.H.  Cohen,  T.  Williams,  K.  Sigvardt,  A.  McClennan, 
S.  Laederich,  S.  Alfords,  M.  Baker,  N.  Mellon,  M.  Dobrof,  S.  Strogatz,  E. 
Marder,  S.  Giszter. 
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COUPLED  OSCILLATORS  AND  THE  DESIGN  OF  CENTRAL  PATTERN  GENERATORS 

N.  Kopell,  Boston  University 

G.B.  Eraen trout,  University  of  Pittsburgh 


Abstract: 

Much  can  be  deduced  about  the  behavior  of  chains  of  oscillators  under 
ainiaal  assuaptions  about  the  nature  of  the  oscillators  or  the  coupling.  This 
paper  reviews  work  on  such  chains,  and  provides  a  fraaework  within  which 
implications  say  be  drawn  about  the  neural  networks  that  govern  undulatory 
locoaotion  in  lower  vertebrates. 
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1.  INTRODUCTION 


One  strategy  in  the  study  of  neural  networks  is  to  attempt  to  find  every 
cell  which  parti pa tea  in  the  network,  and  learn  all  of  the  major  properties  of 
such  cells  and  their  interactions.  Such  a  strategy  is  widely  espoused  [1]. 

For  small  networks,  as  occur  in  some  invertebrates  [2-4],  this  is  a  feasible 
program,  which  is  meeting  some  success.  For  large  networks,  such  as  those 
governing  even  relatively  simple  processes  in  vertebrates,  it  is  unlikely  that 
that  it  can  be  carried  out  any  time  in  the  foreseeable  future;  another 
strategy  is  required.  This  is  particularly  true  when  the  network  in  question 
consists  of  elements  that  have  complex  dynamics  such  as  oscillations,  so  that 
the  behavior  of  the  ensemble  is  not  easily  predictable  from  knowledge  of  the 
components  and  their  connections. 

Fortunately,  much  of  the  behavior  of  nonlinear  oscillators,  when  forced 
or  coupled  to  one  another,  is  independent  of  at  least  some  of  the  detail  of 
the  structure  of  the  oscillators  and  their  coupling.  Indeed,  mathematics 
provides  some  tools,  and  a  framework,  for  sorting  out  which  details  are  in 
fact  important;  within  this  framework,  it  is  possible  to  draw  implications 
from  the  kind  of  phenoamno logical  observations  that  are,  or  potentially  may 
be,  available. 

One  area  in  which  this  point  of  view  can  be  useful  is  the  study  of 
central  pattern  generators  (CPQ’s) .  These  are  neural  networks  thought  to 
govern  a  range  of  stereotypic  motions,  including  such  rhythmic  activities  as 
walking,  running,  swiaiing,  chewing  and  breathing.  It  is  widely  believed  that 
they  contain  neurons  or  neural  subnetworks  that,  in  the  appropriate 


physiological  conditions,  are  capable  of  endogenous  oscillations  [5,6].  The 
mathematical  question  to  be  understood  is  how  these  oscillators  self~organize 
to  fora  the  appropriate  spatio-teaporal  pattern  of  signals  to  the  muscles. 

We  shall  concentrate  on  the  CPG’s  of  the  lamprey  and  other  fish-like 
animals.  The  lamprey  CPQ  is  (me  of  the  most  studied,  since  it  is  a  vertebrate 
having  a  relatively  simple  spinal  cord  which  can  be  dissected  out  and  kept 
alive  for  further  study  [7].  It  provides  a  preparation  in  which  rhythmic 
output  can  be  elicited  even  when  all  influence  of  sensory  feedback  and 
descending  control  from  the  brain  are  removed.  Furthermore,  the  relatively 
simple  geometry  of  the  spinal  cord  also  helps  in  the  analysis. 


2.  THE  FISH  SWPMINQ  CPQ. 

Fish  of  many  species  propel  themselves  through  water  by  rhythmic 
undulations  [8].  The  undulations  are  caused  by  contractions  that  pass  down 
the  axial  muscles  (with  muscles  on  the  opposite  sides  of  the  fish  180°  out  of 
phase).  These  contractions  are  in  turn  directed  by  activity  in  motoneurons, 
whose  outputs  emerge  from  special  positions  along  the  spinal  cord  called 
"ventral  roots"  that  are  arranged  periodically  along  the  cord.  Measurements 
of  the  activity  from  these  ventral  roots,  and  from  innervated  muscles,  have 
been  made  in  several  species.  They  have  been  found  to  be  qualitatively  similar 
to  each  other  and  to  the  measurements  from  the  isolated  spinal  cord  of  the 
lamprey  [8,9]. 

The  measurements  taken  at  a  ventral  root  show  rhythmic  voltage  changes 


(bursts  of  activity)  at  a  uniform  frequency,  with  a  phase  lag  averaging  IX  per 
segment  between  any  two  roots;  the  lags  are  approximately  proportional  to  the 
distance  between  the  roots.  Thus,  there  are  travelling  waves  of  electrical 
activity  having  an  approximately  constant  wave  speed.  (Recent  data  involving 
finer  measurements  suggest  that,  at  least  in  the  isolated  cord,  the  phase  lags 
per  segment  are  not  so  constant  along  the  cord  [10,11]).  The  swimming  speed 
of  an  intact  fish,  which  is  proportional  to  the  frequency  of  oscillation,  may 
vary  over  an  order  of  magnitude  [8].  A  very  striking  observation  is  that  the 
wavelength  of  the  electrical  and  resulting  mechanical  waves  remains  constant, 
at  about  one  body  length,  independent  of  the  speed  of  swimming,  and  this  is 
true  as  well  of  the  isolated  spinal  cord  [9] .  Though  this  is  obviously  useful 
for  the  mechanical  and  hydrodynamical  properties  of  the  animal,  it  is  not 
a-priori  clear  how  such  regulation  is  accomplished,  and  it  is  one  of  the 
achievements  of  the  mathematics  to  suggest  how  this  can  be  done. 

The  work  described  in  this  paper  is  addressed  to  the  question  of  how  the 
network  self-organizes  in  order  to  help  produce  the  appropriate  phase-lags. 
(The  lack  of  complete  constancy  of  the  phase  lags  in  the  isolated  cord 
suggests  that  the  mechanics  and  sensory  feedback  may  also  be  important  in 
regulating  the  lags;  this  requires  further  work.)  The  above  observation  about 
constancy  of  the  wave-length  constrains  the  possible  models.  For  example,  two 
mechanisms  that  are  ruled  out  have  waves  of  activity  due  to  transmission  delay 
or  synaptic  delay.  Since  those  delays  are  generally  time  independent,  if  the 
frequency  of  the  oscillations  is  increased,  the  delay  will  occupy  a  larger 
fraction  of  the  oscillator  cycle,  and  the  phase  lags  between  any  two  points 
will  increase;  this  has  the  effect  of  decreasing  the  wave-length.  We  also 
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note  that  the  constancy  of  the  wave  length  is  a  unusual  property  of  nonlinear 
waves;  it  is  absent  in  the  waves  of  reaction-diffusion  equations  generated  by 
local  oscillations  [12]  or  waves  in  nerve-conduction  equations  [13].  Another 
constraining  observation  is  that  the  wave-producing  aechanisa  is  highly 
distributed;  any  sufficiently  long  section  of  the  spinal  cord  can  "fictively 
swim”.  Finally,  an  intact  animal  is  capable  of  temporarily  swiasing  backwards 
(in  which  case  the  electrical  wave  passes  from  tail  to  head),  and  this  is 
occasionally  seen  even  in  an  isolated  spinal  cord  [9,11]. 

The  cells  participating  in  the  CFG  have  been  only  partially  identified 
[14,15].  Nevertheless,  there  is  such  indirect  evidence  for  the  existence  of 
coupling  in  both  the  rostral  (toward  the  head)  and  caudal  (toward  the  tail) 
directions.  The  aoat  compelling  evidence  involves  the  recent  experiments  of 
Sigvardt  and  McClellan  [16],  who  exploited  the  existence  of  mechanoreceptors 
in  the  spinal  cord  to  force  the  electrical  activity  at  either  the  rostral  or 
caudal  end.  By  cutting  the  direct  connections  from  the  mechanoreceptor  to  all 
but  the  very  end  that  was  being  mechanically  forced,  they  insured  that  any 
longer  range  effects  were  through  the  intersegmental  coordinating  system;  they 
found  that  forcing  at  either  end  of  a  long  portion  of  the  cord  could  entrain 
the  oscillations  of  the  entire  piece,  thus  proving  that  the  coordinating 
system  provides  signals  in  both  directions.  In  other  evidence,  work  of  N. 

Dale  shows  that  if  only  a  portion  of  the  cord  is  chemically  excited,  the 
excitation  can  pass  either  caudal  or  rostral  [14].  Finally,  when  a  portion  of 
a  cord  is  cut  into  two  pieces,  the  frequency  of  the  pieces  are  in  general 
different  from  that  of  the  larger  piece,  suggesting  that  that  the  removed 
connections  had  effects  in  both  the  rostral  and  caudal  directions  [10]. 


For  the  above  reasons,  any  mathematical  description  of  the  CPG  oust  take 
into  account  the  two-way  coupling.  The  laaprey  cord  is  known  to  have  several 
coupling  systeas,  operating  at  different  distances  [17].  We  shall  concentrate 
on  the  short  range  coupling,  partly  because  we  believe  that  it  will  turn  out 
to  be  the  aost  iaportant  for  the  self-organization,  but  also  because  a 
thorough  investigation  of  the  properties  and  limitations  of  short-range 
coupling  can  give  clues  to  the  reasons  for  the  existence  of  the  longer  range 
systeas.  Further  hypotheses  on  the  coupling  will  be  discussed  in  aore  detail 
in  Section  4. 

We  shall  treat  the  network  as  a  chain  of  discrete  oscillators.  (A  aore 
coaplete  description  for  soae  animals,  including  laaprey,  would  be  a  pair  of 
such  chains.  However,  since  the  analogous  oscillators  on  opposite  sides  of 
the  cord  generally  operate  in  antiphase,  it  usually  suffices  to  consider  just 
one  chain. )  We  note  here,  however,  that  the  existence  of  discrete 
oscillators,  each  presumably  (though  not  necessarily)  associated  with  a 
segment  in  the  spinal  cord,  has  not  been  established  in  the  laaprey.  In  the 
leech,  the  spinal  cord  is  segmented,  and  a  single  segment  has  been  shown  to  be 
capable  of  oscillation  [18];  the  lamprey  spinal  cord  is  not  itself  segmented, 
and  no  piece  smaller  than  3  segments  (of  the  associated  musculature)  has  yet 
been  shown  to  be  competent  to  oscillate.  Extensions  of  the  work  on  discrete 
oscillators,  discussed  in  Section  50,  go  in  the  direction  of  showing  how  an 
unsegmented  array  of  oscillators  can  behave  like  a  discrete  chain. 


3.  A  PAIR  OF  SIMPLE  MATHEMATICAL  DESCRIPTIONS. 


A.  Sine-coupled  phase  oscillators. 

The  earliest  attempt  to  aodel  some  of  the  features  of  the  fish  swimming 
C.P.Q.  was  based  on  the  idea  that  a  gradient  in  the  natural  frequencies  of  the 
oscillators  of  a  chain  can  cause  the  formation  of  travelling  waves  in  the 
chain  [19].  To  show  this,  Cohen  et  al.  used  a  simple  and  completely  solvable 
set  of  equations.  The  oscillators  are  represented  by  phases  0^,  k  = 

1....N+1,  which,  in  the  absence  of  interactions,  satisfy  the  equations  9£  = 
u^,  where  is  the  frequency  of  the  oscillator.  The  full  coupled 

equations  are 


•k =  “k  ♦ B  si"<‘k.i  -  V  *  >  *in(Vi  •  V 


where  B  is  a  measure  of  the  strength  of  coupling.  Travelling  waves  are 
phaselocked  solutions  to  (1),  i.e. ,  ones  for  which  the  phase  differences  = 
0k+l  -  0^  are  independent  of  time;  for  systems  of  the  form  (1),  this  is 
equivalent  to  solutions  for  which  the  frequency  9£  is  the  same  for  each  k 
and  is  independent  of  time.  The  equations  for  the  phase  differences  in 

a  phaselocked  solution  then  satisfy 


0  =  d  +  B  M  sin  T 


where  £  =  (^,...fN),  d_=  (dj, . . .  ,djj) ,  d^  is  the  frequency  difference 

-  w^,  sin  f  =  (sin  P^,...,sin  y  )  and  is  the  N  x  N  tri-digonal  matrix 
with  -2  along  the  main  diagonal  and  1  along  the  other  two  diagonals. 


Since  N  is  invertible,  (2)  can  be  solved  for  any  choice  of  {d.},  provided 
/■**  1 

that  every  component  of  B  is  less  than  one  in  absolute  value.  It 

can  be  shown  that  this  requireaent  iaplies,  for  a  gradient  (in  which 
all  have  the  saae  sign),  that  the  aust  scale  like  Od/N2),  i . e . ,  the 

total  frequency  difference  that  can  be  sustained  and  still  have  phaselocking 
scales  like  0(1/N). 

It  is  easy  to  show  that,  in  the  presence  of  a  gradient,  the  solutions  to 
(2),  if  they  exist,  do  not  correspond  to  constant  speed  waves,  i.e.,  solutions 
with  Pj'  independent  of  k  and  non-zero  [19].  If  the  oscillators  are 
identical,  the  solution  is  synchrony,  i.e.,  a  0;  thus  frequency 
differences  are  required  to  get  any  travelling  wave  at  all. 

B.  Chains  with  one-wav  coupling 

Another  way  to  couple  oscillators  in  a  chain  so  as  to  get  a  systes  that 
is  easy  to  analyze  is  to  aasuae  that  the  coupling  is  only  in  one  direction. 

For  such  a  chain,  the  analysis  of  the  enseable  reduces  to  the  analysis  of  any 
two  successive  oscillators,  with  the  first  driving  the  second,  which  drives 
the  third,  etc.  If  the  forcing  is  done  by  a  localized  stimulus,  such  as  a 
rapid  pulse,  and  if  the  resulting  trajectories  return  very  quickly  to  the 
liait  cycle,  the  behavior  of  each  pair  of  forced  oscillators  can  be  analyzed 
using  "phase  response  curve"  theory  [20,21]. 

A  phase  response  curve  for  a  forced  oscillator  gives  the  change  in  phase 
(after  the  transients  have  died  away)  due  to  the  stiaulus  as  a  function  of  the 
phase  at  which  the  stimulus  is  given.  Fron  such  a  curve,  it  is  possible  to 
calculate  at  what  frequencies  the  stimulus  say  be  given  in  order  to  entrain 


the  oscillator.  In  general  the  phase  of  the  forced  oscillator  at  which  the 
stimulus  appears  when  there  is  entrainaent  varies  with  the  frequency  of  the 
forcing,  and  can  also  be  calculated  fros  the  PRC  [22].  Even  if  the  forcing 
oscillator  is  identical  to  the  forced  one,  the  phase  difference  between  the* 
need  not  be  zero. 

A  chain  of  oscillators,  each  forcing  its  successor  by  a  single  pulse  and 
whose  transients  are  very  rapid  aa  above,  can  be  analyzed  using  PRC  theory. 

One  sees  iaaediately  that  such  a  chain  can  support  a  constant  speed  (i.e. , 
constant  phase-difference)  travelling  wave,  provided  that  the  oscillators  are 
identical.  Furthermore,  the  wave  can  go  in  either  direction,  i.e.,  it  needn’t 
be  in  the  saae  direction  as  the  coupling;  the  direction  depends  on  the  shape 
of  the  PRC.  This  direction  is  unchanged  by  the  presence  of  a  snail  frequency 
gradient  which,  for  sine-coupling  discussed  above,  would  sake  the  waves  travel 
in  the  opposite  direction.  However,  for  a  sufficiently  large  adverse 
gradient,  the  wave  can  reverse  direction.  (For  a  still  larger  gradient,  there 
is  no  locking  at  all.)  Such  a  aodel  has  been  applied  to  understand  the 
so-called  "aetachronal  waves"  in  the  swiaaerets  of  crayfish  [23]. 

Both  of  the  above  descriptions  are  useful  for  understanding  some  of  the 
properties  of  the  fish  swiaaing  CPG,  but  neither  is  wholly  adequate.  As  was 
discussed  in  Section  2,  the  coupling  is  in  both  directions,  not  a  unilateral 
forcing.  Although  the  PRC  fraaework  can  be  used  to  describe  two-way  coupling, 
the  ease  of  analysis  vanishes,  and  solutions  aust  be  found  nuaerically.  More 
crucially,  the  PRC  fraaework  deals  with  forcing  or  coupling  by  pulses;  it 
cannot  deal  with  coupling  that  is  aore  spread  out  within  each  cycle,  with 


■ultiple  pulses  per  cycle  [4]  or  graded  release  of  transaitters  [24,25]. 
Indeed,  there  is  Mathematical  evidence  that,  unless  the  coupling  is  spread  out 
as  above,  then  if  the  coupling  if  two-way,  there  nay  be  no  rhythaic  activity 
at  all.  This  is  discussed  in  More  detail  in  Section  4F  and  5D. 

The  aajor  difficulty  with  the  first  description  is  the  specificity  of  the 
coupling;  it  is  difficult  to  tell  a-priori  which  of  the  conclusions  that  are 
drawn  depend  on  this  particular  coupling.  As  will  be  seen  below,  the  features 
that  aake  this  exaaple  easy  to  analyze  also  aake  it  sufficiently  "non-generic" 
that  it  is  a  bed  guide  to  intuition.  (For  exanple,  the  conclusion  that  locking 
fails  for  large  N  unless  the  total  frequency  differences  scale  like  1/N 
turns  out  to  be  generically  wrong.  Other,  nore  qualitative  differences  are 
discussed  in  Section  4B.)  We  shall  present  below  a  framework  that  includes 
the  above  examples,  while  avoiding  the  undesirable  restrictions. 

4.  A  VERY  GENERAL  MATHEMATICAL  FRAMEWORK 

At  TIW 

Since  the  composition  of  any  fish  swiaiing  network  is  largely  unknown, 
and  the  networks  can  be  expected  to  vary  aaong  species,  it  is  desirable  to 
build  a  theory  that  is  as  robust  as  possible.  For  the  oscillators  themselves , 
our  framework  is  alaost  completely  general.  Bach  oscillator  can  be  any 
system,  a  single  neuron  or  a  subnetwork,  that  can  be  described  by  an  ordinary 
differential  equation  having  a  stable  periodic  limit  cycle.  The  state 
variables  of  these  equations  can  be  vectors  of  arbitrary  size,  including  any 


variables,  such  aa  membrane  potentials,  firing  rates,  channel  openings,  etc., 
that  may  be  relevant.  The  oscillators  at  different  positions  along  the  chain 
need  not  be  the  setae  (for  example,  they  need  not  have  the  same  frequency),  but 
the  changesnust  be  gradual  along  the  chain.  The  large  aaount  of  generality 
about  the  oscillators  reflects  the  fact  that  the  behavior  of  the  system  is 
insensitive  to  the  details  of  the  oscillators. 

By  contrast,  as  will  be  seen  in  this  paper,  soae  of  the  details  of  the 
coupling  are  quite  important.  We  shall  start  with  a  few  hypotheses;  these  are 
needed  to  carry  out  the  analysis,  but  we  also  believe  that  they  represent 
important  constraints  on  the  design  of  the  network,  as  we  argue  below  in 
Section  4F  and  5D.  We  shall  not  assume  that  the  coupling  among  the 
oscillators  is  only  between  nearest  neighbors;  indeed,  the  coupling  to 
somewhat  more  distant  neighbors  turns  out  to  have  important  design 
implications.  (See  Sections  40  and  5A. ) 

We  write  below  the  form  of  the  equations  to  be  discussed.  For 
simplicity,  they  are  written  only  for  nearest  neighbor  coupling;  if  there  is 
additional  coupling,  more  terms  must  be  added. 

(3)  =  W  *  *kQ*(Xk+l,Xk)  +  Vl^VrV’ 

In  (3),  X£  =  F^(X^)  represents  the  equation  of  the  oscillator  at  the  k*'*1 
position.  The  {s^}  are  positive  numbers  representing  the  strength  of  the 
coupling,  which  is  allowed  to  vary  (gradually)  along  the  chain.  s^G*  and 
*k-lQ  represent ,  respectively,  ascending  and  descending  coupling,  which  need 
not  be  the  same.  To  add  coupling  with  other  neighbors,  we  include  terms  of 
.Xj'),  multiplied  by  appropriate  strengths. 


the  torn 


So  far  the  coupling  ia  almost  completely  general.  We  now  add  two 
hypotheses.  The  first  one  is  that  all  the  state  variables  of  the  oscillators 
are  completely  determined  by  the  phases  of  the  oscillators.  This  ia  always 

true,  provided  that  the  coupling  is  not  so  strong  as  to  overcome  the  tendency 

of  limit  cycle  oscillators  to  maintain  their  wave-form  and  amplitudes,  and  is 
a  fairly  benign  hypothesis.  (See  (26}  for  more  details.) 

The  second  is  by  no  means  so  automatic,  and  is  therefore  much  more 
interesting.  It  is  assumed  that  if  the  coupling  terms  in  (3)  are  replaced  by 

their  averages  over  a  cycle  of  the  k**1  oscillation,  then  the  resulting 

equations  have  the  same  phaselocking  behavior  as  (3).  If  the  coupling  is 
"weak",  it  follows  from  so-called  "averaging  theory”  that  this  hypothesis  can 
be  made  without  losa  of  generality  [27];  for  large  coupling,  it  need  not  hold. 
However,  if  the  coupling  signals  are  sufficiently  dispersed  around  the  cycle, 
instead  of  being  concentrated  in  a  small  "window”  of  phases,  the  hypothesis 
can  hold  even  for  strong  coupling  (28).  Furthermore,  as  mentioned  in  Section 
3B,  current  numerical  and  analytical  evidence  indicates  that,  when  coupled  in 
both  directions  by  pulses,  the  chain  (or  even  a  pair)  of  oscillators  may  stop 
all  rhythmic  activity.  [28].  Thus,  we  believe  that  this  hypothesis  is  more 
than  a  mathematical  convenience.  Other  hypotheses  to  be  discussed  below  also 
emerge  as  conclusions  of  the  mathematical  investigation,  i.e. ,  they  are 
further  restrictions  that  are  necessary,  within  this  framework,  to  get  the 
observed  behavior,  and  therefore  also  constitute  possible  restrictions  on 
design. 

* 

Under  the  above  hypotheses,  the  full  equations  (3)  can  be  reduced  to  a 
much  smaller  set  of  equations  describing  the  interactions  of  the  phases  of  the 


hi 


oscillators  [26,29].  If  0^  is  the  phase  of  the  k**1  oscillator,  then  after 
averaging,  these  equations  have  the  fora  : 

(4)  tfk  =  *  *kH  (*k+l  “  V  *  Sk-1H  (*k-l  '  V 

[26,29].  Here,  H+  and  H  are  scalar,  2 n- periodic  functions,  and  is  the 
frequency  of  the  k^  oscillator.  The  ability  to  reduce  to  phase  equations 
is  the  first  hypothesis;  the  fact  that  these  equations  have  the  fora  (4),  with 
the  coupling  teras  dependent  only  on  the  differences  of  the  phases,  is  a 
consequence  of  the  averaging  procedure.  Note  that  equations  (2)  are  of  this 
fora,  with  s^  a  1,  H+(*)  =  H  (•)  =  sin(*).  We  shall  discuss  below  the 
behavior  of  the  nearest  neighbor  equations  (4),  and  return  in  Section  G  to  the 
■ore  general  case. 


otonic  v 


No  further  hypotheses  are  needed  in  order  to  be  able  to  analyse  the  above 
equations.  However,  because  of  the  generality  of  the  fraaework,  it  becosres 
possible  to  see  the  importance  of  certain  further  distinctions.  We  focus  now 


on  the  aatheaatical  differences  between  electrotonic  coupling  and  coupling  due 
to  chemical  synapses. 


For  electrotonic  coupling  (or,  equivalently,  bulk  diffusion  coupling), 
the  coupling  teras  G±(X,X)  vanish  if  X  =  X;  for  the  averaged  equations, 
this  implies  that  H*(0)  =  0.  By  contrast,  if  two  neur,onal  oscillators  are 
coupled  via  the  release  of  transaitter  substances,  they  can  continue  to 
influence  one  another  even  when  they  are  in  the  saae  state,  and  they  can,  on 


the  average,  slow  each  other  down  or  speed  each  other  up.  For  the  averaged 
equations,  this  corresponds,  respectively,  to  H^O)  <  0  and  H^O)  >  0. 

(Sven  if  H+(0)  and  H  (0)  have  opposite  signs,  the  frequency  in  a  chain  of  more 


than  two  identical  oscillators  will  not  in  general  be  the  sane  as  the  natural 
frequency.  See  Section  4C. )  We  shall  refer  to  coupling  which  has  the 
property  H*(0)  SO  as  "synaptic  coupling",  with  the  understanding  that  this 
is  a  property  that  can  be  held  by  oscillating  subnetworks  as  well  as 
individual  neurons. 

Synaptic  coupling  is  i sport ant  when  oscillators  are  arranged  in  a  linear 
geoswtry  because  it  forces  the  existence  of  non-zero  phase  differences  aaong 
the  oscillators.  To  see  this,  we  let  a  ~  *jl»  and  use  (4)  to  get 


<s*>  ViH  <Vi>  -  *kB  <V  *  -  ViH~(“Vi> 


j  -  the  frequency  difference  between  successive  oscillators. 
Because  the  oscillators  on  the  ends  are  connected  to  only  one  neighbor  each, 
the  first  and  last  equations  of  (5a)  differ  fro«  the  above  in  sissing  the 
first  or  last  tens;  since  the  {s^}  are  positive,  this  eeans  that 


Here 


\  = 


(5b) 


h-(-V  ■ 0  = 


Phase locked  solutions  to  (3)  or  (4)  correspond  to  critical  points  of  (5), 
i.e. ,  solutions  in  which  the  phase  differences  are  independent  of  tine.  If 
there  is  diffusive  coupling  of  identical  oscillators  (d^  =  0  for  all  k),  and 
the  strengths  of  coupling  do  not  depend  upon  k,  the  network  oscillates 


synchronously, i. e. ,  =  0  for  all  k.  (Equation  (2)  can  be  derived  froa 

equations  of  the  fora  (3)  having  diffusive  coupling.)  When  the  coupling  is 
synaptic,  synchrony  is  iapossible  even  for  the  above  sinple  case  4^  s  0,  s^  s 
1:  it  is  easily  seen  that,  because  of  (5b),  =  0  does  not  solve  the  first 

or  last  equation  of  (5a).  Thus,  synaptic  coupling  forces  non-zero  boundary 
conditions.  As  will  be  seen  below,  this  effects  the  behavior  far  away  froa 
the  edges.  We  believe  that  synaptic  coupling  is  the  priaary  aechanisa  for  the 
production  of  the  waves  of  electrical  activity,  as  we  explain  further  in 


Section  5. 


We  reaark  that  the  synaptic  coupling  condition  is  the  analogue  for  (4)  or 
(5)  of  the  aechanisa  in  Section  3B  that  produced  produced  travelling  waves  (p^ 
P  0)  even  with  identical  oscillators.  That  is,  consider  a  pair  of  oscillators 
having  period  1,  coupled  in  one  direction  by  pulses,  and  let  6=0  denote 
the  phase  at  which  the  coupling  pulse  is  eaitted.  Then  a  non-zero  phase 
difference  p  is  produced  if  the  PRC  function  4(6)  has  a  zero  6fl  P  0 
with  0  <  d'(60)  <  2  for  stability;  if  4(0)  =  0  and  0  <  4‘(0)  <  2, 
synchrony  will  occur.  If  the  coupling  by  pulses  is  not  too  strong,  then 
averaging  theory  nay  be  applied,  and  a  function  H  coaputed  as  in  (4).  It 
can  then  be  shown  that  the  condition  H  (0)  P  0  is  equivalent  to  .4(0)  *  0. 
(See  [22]  and  [30]  for  aore  details.) 


It  is  difficult  to  find  analytically  even  the  tiae- independent  solutions 


to  (5).  However,  it  can  be  shown  [29,31]  that  the  latter  solutions  are 
arbitrarily  well  approximated,  for  N  large,  by  solutions  to  an  associated 
continuua  equation 


(6a) 


0  =  /J(x)  +  (a(x)f(P)]x  +  ^[s(x)g(^))xx. 

Here,  0  <  x  <  1,  and  ?(x)  is  a  continous  function  whose  discretization  is 
the  unknown  phase  differences  i*®**  ^  **  *(k/N).  Siailarly,  if  = 

\/V,  then 

/S(x)  and  s(x).  The  2 w- periodic  scalar  valued  functions  f  and  g  are 
computed  froa  H+  and  H  ;  for  example,  if  H+  =  H  a  H,  then  f  is  the  even 
part  of  H  and  g  is  the  odd  part  of  H.  For  a  relatively  simple  class  of 
oscillators  and  such  isotropic  coupling,  f  and  g  have  the  form  f  = 

A  cos(P)  and  g  -  B  sin(r),  for  same  A,B  [26].  In  the  anisotropic  case,  f 
and  g  are  related  to  H+  and  H  by  (f+g)(*)  =  H+(P)  and  (f-g)(P)  = 

H  (~r) .  (See  [29]  for  details.)  The  end  conditions  (5b)  translate  into 
boundary  conditions 


r  *  rL  at 

x  =  0, 

where 

H~(-rL)  =  0 

X  =  1, 

where 

hVr)  -  0. 

Equations  (6)  fora  a  two-point,  singularly  perturbed  boundary  value  problea. 
Such  problems  need  not  have  any  solutions,  and  solutions  need  not  be  unique. 
Nevertheless,  it  can  be  shown  for  these  equations  that  if  the  range  of  natural 
frequencies  {w^}  is  not  too  large,  and  some  not  very  restrictive  technical 
hypotheses  are  satisfied,  then  there  is  indeed  exactly  one  solution,  and  there 
is  a  stable  solution  to  (5)  which  converges  to  this  as  N  grows  without 
bound.  One  of  these  hypotheses  is  a  "genuine  nonlinearity"  condition,  i.e., 
f"  s  0  over  the  part  of  its  domain  where  the  solution  turns  out  to  lie. 


and 


<\> 


are  discretizations  of  s 


continuous  functions 


The 


sign  of  f"  plays  a  role  discussed  below.  We  note  that  the  continuum 
equation  (6)  is  not  obtained  froa  (5)  by  a  standard  continuua  limit,  but  is 
■ore  closely  related  to  the  idea  in  physics  of  a  "thermodynamic  limit".  (See 
(29]  for  more  details. ) 

The  importance  of  the  equations  (6)  is  that  they  are  much  easier  to  solve 
analytically  than  (5),  and  it  is  even  easier  (though  often  still  not  trivial) 
to  get  approximate  solutions  to  (6).  Thus  it  is  possible  to  use  these 
equations  as  a  kind  of  "calculus"  to  figure  out  the  properties  of  the 
solutions  to  (5).  To  do  this,  we  note  that  many  singular  perturbation 
problems,  including  (6),  have  solutions  that  behave  over  much  of  the  domain  as 
if  the  small  term  (in  this  case  (1/N) (s(x)g(P) ]  )  were  absent.  Since  the 

omission  of  that  term  leaves  a  first  order  equation  that  cannot  solve  both 
boundary  conditions,  there  is  a  region,  called  a  transition  layer,  in  which 
that  term  is  significant.  Such  a  layer  is  usually,  but  not  always,  near  a 
boundary,  and  is  then  called  a  "boundary  layer". 

If  the  frequencies  are  monotone  (so  /8(x)  >0  or  $  0  for  all  x)  and  the 
{Sj'}  are  constant,  there  is  generically  a  boundary  layer  at  only  one  end  [29]. 
Then  the  solution  is  very  close  over  most  of  the  domain  to  the  solution  to  the 
"outer  equation" 

(7))  0  =  /}(x)  +  (s(x)f(r)]x 

satisfying  the  boundary  condition  at  the  other  end.  Under  some  robust 
conditions,  discussed  in  Section  40,  there  may  be  boundary  layers  at  both 
ends;  in  this  case  the  solution  to  (6)  is  very  close  to  a  solution  to  (7) 
satisfying  a  condition  at  an  interior  point  instead  of  at  a  boundary  [31]. 


The  choice  of  boundary  condition  to  be  satisfied  by  the  outer  equation 
determines  much  of  the  solution.  For  example,  if  the  oscillators  all  have  the 


same  natural  frequency,  and  the  strengths  s^  are  independent  of  k,  then  the 
solutions  r(x)  to  (7)  are  constant  functions.  Thus  the  boundary  condition 
determines  the  value  of  the  constant  function.  If  that  boundary  value  is 
negative,  so  =  0^.^  -  0^  <  0,  there  is  a  wave  that  goes  in  the  direction 
of  increasing  k  (forward  swimming);  if  the  boundary  value  is  positive,  the 
wave  goes  in  the  direction  of  decreasing  k  (backward  swimming).  Furthermore, 
if  the  boundary  conditions  are  of  opposite  signs,  as  is  true  if  H+  and  H 
are  multiples  of  one  another,  then  the  position  of  the  boundary  layer 
determines  the  direction  of  swimming.  (The  boundary  conditions  need  not  have 
opposite  signs,  in  which  case  the  boundary  layer  may  be  much  smaller,  or  even 
nonexistent.)  We  also  note  that  if  max|H*|  0,  so  that  the  coupling  becomes 
only  descending,  then  the  outer  solutioo  approaches  the  full  solution  to  the 
uni-directional  coupling  problem  [31]. 

The  boundary  condition  to  be  satisfied  also  determines  the  frequency. 

For  example,  if  (6b)  is  satisfied  at  x  =  0,  it  can  be  shown  that  the 


phase locked  frequency  converges,  as  N  -»  <*,  to  the  right  hand  side  of  (4), 


with  k  =  1  and  0^  -  0^  replaced  by  a  similar  computation  holds  for 

the  other  cases.  Note  that,  in  this  case,  the  sign  of  H+(f^)  determines 
whether  the  frequency  of  the  ensemble  will  be  higher  or  lower  than  the 
frequency  of  the  first  oscillator;  when  the  coupling  is  synaptic,  it  will  not 
be  the  same  as  the  latter.  (If  the  boundary  condition  at  x  =  1  is  satisfied, 

it  is  the  sign  of  H  ( - that  is  relevant.) 

When  there  is  only  one  boundary  layer,  its  position  is  determined  in  a 


V  V 


complicated  way  from  the  signs  of 

A 


f  and  Q  =  2[f(fL)  - 


f (i»_)  ]  -  /  /}(y)dy,  where  the  latter  measures  a  balance  between  the  overall 
frequency  gradient  and  the  amount  of  anisotropy.  If  H+  and  H  are 
multiples  of  one  another,  and  there  is  no  frequency  gradient,  it  can  be  shown 
that  for  f • •  >0,  it  requires  stronger  coupling  in  the  forward  direction  to 
have  forward  moving  waves;  for  f ' '  <0,  the  opposite  is  true.  (See  [29].) 


D.  The  effect  of  "pacemakers"  and  "lesions" 

If  coupling  strengths  are  constant  and  frequencies  are  monotonic,  then 
the  outer  solution  (almost)  always  satisfies  one  of  the  two  boundary 
conditions.  (The  exceptions  are  special  cases,  such  as  identical  oscillators 
and  H+  =  H  ,  which  contain  symmetries  that  can  be  perturbed  away.)  However, 
if  there  is  a  region  in  which  the  frequencies  are  sufficiently  different  from 
the  ambient  local  frequencies  a  new  phenomenon  may  occur.  For  simplicity,  we 
describe  the  almost-isotropic  case;  the  generalizations  are  in  [31],  If  the 
local  frequency  is  sufficiently  higher  (if  f"  >  0,  or  lower  if  f"  <  0),  then 
the  outer  solution  satisfies  a  condition  at  the  point  of  maximum  frequency 
instead  of  at  one  of  the  boundaries,  and  there  are  in  general  boundary  layers 
at  both  ends.  (See  Figure  la-c).  In  this  case,  waves  propagate  outwards  from 
the  point  of  highest  frequency  (if  f"  >  0,  or  inwards  if  f”  <  0.).  In  (1), 
corresponding  to  diffusive  coupling,  this  "pacemaker  effect"  occurs  whenever  a 
state  of  uniform  frequency  is  perturbed  by  local  increase.  When  there  is 
synaptic  coupling,  in  general  there  is  a  threshold  for  this  effect  to  happen; 
for  smaller  local  increases  in  local  frequency,  the  direction  of  the  wave  is 
unchanged,  so  the  synaptic  coupling  acts  to  buffer  the  wave  direction  against 


perturbations  in  frequency.  The  frequency  of  the  ensemble  is  also  unchanged 
(to  lowest  order  in  1/N)  until  the  threshold  is  reached,  and  then  it  changes 
dissent iauemsiy  (in  the  limit  as  N  -»  «).  Thus,  the  crossing  of  this 
threshold  corresponds  to  a  kind  of  phase  transition  for  the  solution  [31]. 

For  sufficiently  large  local  increases  or  decreases  in  frequency,  locking  is 
lost. 

Similar  effects  are  found  if  instead  of,  or  in  addition  to,  nonmonotonic 
frequencies,  we  allow  the  coupling  strengths  to  be  locally  reduced.  (This 
corresponds  to  lesions  in  the  spinal  cord. )  However,  the  effect  of  such  a 
"lesion"  depends  greatly  on  the  context;  for  example,  in  a  chain  of 
oscillators  having  a  gradient  in  frequency,  the  same  lesion  could  act  to 
either  locally  increase  or  locally  decrease  the  phase  differences,  depending 
on  the  position  of  the  chain  at  which  the  lesion  is  made.  For  large  enough 
lesions  (but  not  too  large  to  prevent  locking),  the  waves  propagate  outwards 
from  the  lesion  as  for  pacemakers.  See  [31]  for  more  details. 

F.  Averaging  and  oscillator  death 

We  return  now  to  the  hypothesis  introduced  in  Section  3  that  the 
phaselocking  behavior  of  the  chain  of  oscillators  remains  qualitatively  the 
same  if  the  coupling  terms  of  the  k^*1  equation  are  replaced  by  ones  that 
are  averaged  over  the  cycle  of  the  kt^1  oscillator.  As  mentioned  above,  this 
hypothesis  does  not  hold  automatically,  except  when  the  coupling  is 
"sufficiently  weak".  The  natural  question  is  then:  what  is  the  behavior  when 
the  hypothesis  is  violated? 

Numerical  simulations  and  analytic  examples  indicate  that  when  the 


hypothesis  is  badly  violated,  e.g.,  when  the  coupling  is  by  a  single  strong 
pulse  per  oscillator  per  cycle,  the  systen  say  cease  to  oscillate  [28] .  In 
its  simplest  fora,  this  phenomena  occurs  for  a  pair  of  oscillators  coupled 
symmetrically;  it  is  not  exhibited  by  a  forced  oscillator.  A  simple  example 
[32]  is  as  follows:  Let  9^  and  92  he  phases  of  a  pair  of  oscillators,  and 
suppose  that  each  emits  a  pulse  P(9^)  >  0  with  maximal  amplitude  around  0^ 

=  0.  (See  Figure  2.)  The  effect  of  the  this  pulse  is  presumed  to  depend  on 
the  phases  at  which  it  is  received,  that  is,  the  full  equations  are 

(8)  +  a  P(92)  H(9j) 

=  «2  +  a  P(9l)  R(02) 

where  a  is  a  measure  of  the  strength  of  the  coupling.  The  function  R 

plays  the  same  role  as  the  phase  response  curve  described  in  Section  3.  For 

e.g.,  R(9)  =  A  cos(y)  -  B  sin(P),  B  >  0,  and  a  small,  (8)  has  a  stable  limit 

cycle.  It  is  easy  to  show  by  phase  plane  methods  that  for  close  to  one 

another,  as  a  increases,  the  frequency  of  the  limit  cycle  goes  to  zero  and  a 
stable  critical  point  emerges.  (For  =  w2*  ^  *s  evon  ®a®ier  to  check  this 

for  solutions  satisfying  9^  a  6^.) 

This  "death"  of  rhythmicity  appears  to  be  very  robust.  That  is,  it  can 
be  shown  that,  for  large  classes  of  pairs  of  oscillators  whose  interactions 
depend  on  the  phases  as  well  as  the  phase  differences  of  the  oscillators,  a 
sufficient  strength  of  interaction  will  make  the  oscillation  disappear  [28]. 
The  disappearance  involves  the  creation  of  a  stable  time-independent  solution 
for  the  pair,  instead  of  the  stable  limit  cycle.  In  models  that  depend  on 
amplitudes  as  well  as  phases,  such  critical  points  are  also  found.  The  same 
phenomena  occur  in  chains  of  oscillators. 


Since  this  outccae  is  clearly  to  be  avoided  by  networks  that  Bust  compute 
an  oscillatory  pattern,  it  is  natural  to  ask  how  an  oscillatory  network  can  be 
designed  that  will  not  fail  in  the  above  Banner.  The  key  observation  is  that 
such  a  failure  is  due  to  the  non-uniformity  in  local  frequency  along  each  of 
the  cycles;  the  interactions  cause  the  oscillators  to  slow  down  at  sone 
portion  of  the  cycle,  and,  for  large  enough  interactions,  to  stop.  Thus,  a 
design  that  keeps  the  local  frequency  of  each  oscillator  sore  uniform  allows 
for  larger  interactions  without  risking  failure.  Such  unifonsization  can  be 
accomplished  by  having  interactions  between  oscillators  occur  via  coupling 
signals  that  are  dispersed  around  the  cycle,  with  all  (or  at  least  nost)  of 
the  signals  operating  in  a  consistent  way  (i.e. ,  to  phase-advance  or  phase- 
delay  the  receiving  oscillator).  This  is  discussed  further  in  Section  5D. 


Q.  Multiple  coupling 

Moat  of  the  results  described  above  for  nearest  neighbor  coupling  still 
hold  in  the  aore  general  situation  of  aultiple  neighbor  coupling.  However,  we 
believe  that  same  added  effects  of  aultiple  neighbor  coupling  are  of 
considerable  importance. 

Multiple  coupling  is  modelled  by  the  addition  to  (3)  of  terms  of  the  form 
Si(Xk±i*  *k)’  by  appropriate  strengths.  After  averaging,  this 

leads  to  terms  in  (4)  of  the  form  ~  •j<).  As  in  the  nearest  neighbor 

case,  a  continuum  equation  can  be  associated  with  the  discrete  equation  for 
the  phase  differences.  Away  from  the  boundaries,  it  has  been  shown  that  this 
continuum  equation  and  the  discrete  system  behave  in  a  similar  manner  [33]. 

If  K  neighbors  are  connected,  the  continuum  equation  has  the  same  form  as 


£*/ 


(6),  with  t(t)  and  g(r)  replaced  by  f(y)  and  g(r) ,  where 


(g)  f(i»)  =  f  1(r)  +  ...+  fK(K»») 

g(r)  =  «1(r)  +  ...+  gK(R?) 

and  f\,  g^  are  computed  froa  the  connections  to  the  i**1  neighbor.  Under 
very  general  conditions  on  the  {f^}  and  (g^>,  the  resulting  functions  f 
and  g  have  the  property  that  for  the  analogue  of  the  outer  equation  (7),  the 
effect  of  frequency  differences  on  phase  lags  is  considerably  saaller  than 
when  there  is  nearest  neighbor  coupling;  that  is,  the  multiple  coupling  acts 
in  such  a  way  as  to  buffer  the  phase  lags  against  possible  variation  in 
frequency  [33].  It  also  allows  for  larger  frequency  differences  without  loss 
of  coordination  in  the  chain. 

The  boundary  condition  to  be  satisfied  by  the  solution  to  the  new  "outer 
equation"  is  less  understood.  If  (5b)  is  replaced  by  the  analogous 
conditions,  namely,  the  omission  of  all  terms  corresponding  to  oscillators 
with  indices  outside  1,...,N,  then  the  resulting  solution  appears  numerically 
to  satisfy  a  boundary  condition  close  to  that  of  (6b),  with  H+(y)  replaced 
by  (f+g)(y)  and  H  (- t )  replaced  by  (f-g)(V).  (This  results  in  values  ^ 
and  that  are  closer  to  zero  than  for  nearest  neighbor  coupling. ) 

However,  numerical  simulations  and  formal  computations  on  the  discrete 
equations  indicate  that,  when  there  is  multiple  coupling,  boundary  layers  are 
to  be  expected  on  both  edges,  and  these  layers  can  be  oscillatory  [33]. 

In  some  simulations  done  with  a  chain  of  identical  oscillators,  multiply 


coupled,  there  was  a  pronounced  dip  in  the  value  of  the  phase  lag  at  a 
position  corresponding  to  the  nuaber  of  neighbors  coupled.  The  H+  and  H 
for  those  simulations  were  chosen  to  be  compatible  with  the  forcing 
experiments- of  Sigvardt  and  McClellan  [16]  (see  Section  5C).  The  simulation 
results  resemble  the  experimentally  measured  lags,  which  show  such  a  dip  at  a 
distance  of  about  10  segments  from  each  end  of  the  portion  of  the  cord 
measured  [10];  the  effect  appears  to  be  almost  independent  of  the  position  of 
this  portion  within  the  cord.  Both  theory  and  experiment  need  more  work  to 
make  the  results  definitive. 


5.  A  FEW  CONCLUSIONS  AND  SOME  SPECULATIONS 

A.  Phase  lags  and  frequency  differences 

In  Section  4,  we  saw  that  synaptic  coupling  i^>liee  that  there  are 

non-zero  phase  lags,  even  if  the  oscillators  and  the  coupling  strengths  are 

all  the  same.  Indeed,  we  sew  that  the  solution  in  this  cese  has  t .  almost 

k 

constant  (except  near  a  boundary  layer).  This  corresponds  to  constant  speed 
swimming,  with  the  direction  of  the  wave  dependent  on  the  sign  of  the 
Thus,  the  mathematics  shows  that,  if  there  is  synaptic  coupling,  constant 
phase  lags  can  be  achieved  by  regulating  the  frequencies  and  coupling 
strengths  to  be  uniform. 

The  mathematics  also  shows  that  unless  the  frequencies  and  coupling 
strengths  are  uniform,  the  phase  lags  cannot  be  expected  to  be  uniform.  In  a 
very  early  attempt  to  model  the  travelling  waves  in  the  cord,  it  was  suggested 


that  the  waves  were  due  to  a  gradient  in  natural  frequencies,  for  which  there 
was  soae  preliminary  evidence  [19].  According  to  the  above  aatheaatics,  if 
there  is  a  gradient  in  natural  frequencies,  and  if  there  is  not  too  auch 
anisotropy  in  the  coupling,  the  earlier  model  is  correct  in  its  assertion  that 
there  is  a  travelling  wave  starting  froa  the  high  frequency  end.  (This  is 
true  whether  or  not  the  coupling  is  synaptic.  If  there  is  sufficient 
anisotropy,  the  direction  of  the  wave  aay  be  reversed.)  However,  these  waves 
never  have  phase  lags  that  are  independent  of  k,  so  the  existence  of  a 
frequency  gradient  cannot  be  construed  as  the  primary  cause  of  the  constant 
phase  lag  waves.  Indeed,  recent  measurements  of  the  phase  lag  Bade  on 
portions  of  the  lamprey  spinal  cord  show  that  there  are  frequency  differences 
among  different  parts  of  each  isolated  spinal  cord,  with  no  constant  pattern 
in  the  variation  from  one  animal  to  another  [10].  There  is  also  preliminary 
evidence  that  the  phase  lags  vary  less  than  might  be  expected  from  such 
frequency  variation  if  there  were  only  nearest-neighbor  coupling  [10).  These 
results  are  consistent  with  the  conjecture  that  the  lamprey  cord  is  highly 
multiply  coupled,  and  that  the  multiple  coupling  performs  the  buffering 
function  discussed  in  Section  4G.  Simulations  involving  multiple  coupling 
also  predicts  edge  effects,  which  were  seen  in  experimentally  measured  lags 
[10]  (see  section  4G).  We  also  recall  that  synaptic  coupling,  as  opposed  to 
diffusive  coupling,  buffers  the  system  froa  allowing  small  changes  in 
frequency  to  locally  or  globally  change  the  direction  of  wave  propagation. 

For  another  physiological  system  for  which  these  ideas  are  relevent,  see 
[34.]. 


B.  Some  physiological  mechanisms 

It  is  now  easy  to  envision  some  general  physiological  mechanisms  which 
can  accomplish  the  regulation  of  the  phase  lags  as  swimming  speed  (which  is 
proportional  to  the  frequency  of  the  phaselocked  ensemble  [8J )  is  changed.  As 
discussed  above,  the  uniform  phase  lag  that  is  produced  when  the  oscillators 
are  all  the  same  is  determined  by  the  relevant  boundary  condition  (6).  Thus, 
any  change  that  preserves  that  boundary  condition  also  preserves  the  phase 
lags.  In  particular,  if  the  natural  frequencies  are  uniformly  increased, 
without  changing  the  coupling,  there  is  no  change  in  the  equation  or  the 
boundary  conditions  for  the  phase  lags.  Such  uniform  changes  could  be 
produced  by  a  uniform  increase  in  descending  or  ascending  tonic  drive,  and  is 
consistent  with  known  physiology.  We  also  note  that  the  regulation  of  phase 
lags  under  an  increase  in  frequency  for  all  oscillators  can  take  place  even  if 
the  frequencies  are  not  uniform  along  the  cord. 

Similarly,  the  mathematics  suggests  mechanisms  for  changing  the  direction 
of  swumiing.  This  direction  depends  on  the  sign  of  the  If  and 

have  opposite  signs,  this  depends  on  the  position  of  the  boundary  layer,  which 
in  turn  depends  on  the  relative  strength  of  the  forward  and  backward  coupling. 
(See  Section  4.)  Thus,  the  direction  of  swimming  may  be  changed  by  changing 
the  bias  in  the  coupling  that  determines  the  position  of  the  boundary  layer; 
there  is  no  need  for  a  separate  network  to  govern  backward  swimming.  (If  both 
boundary  conditions  and  are  negative,  which  requires  significant 

anisotropy  in  the  coupling,  it  would  take  a  larger  change  to  produce  backward 
swimming;  H  and/or  H  must  be  changed  so  that  and/or  T ^  is 

positive. )  Backward  swiMing  can  also  be  produced  by  appropriate  changes  in 


local  frequency.  The  above  aechanisas  can,  in  principal  be  differentiated 
from  one  another  by  experiaent,  but  the  existing  data  ia  not  sufficiently 
fine.  (See  the  next  section  for  soae  evidence  in  that  direction.) 

C.  Experimental  Manipulations 

The  aatheaatical  fraaework  allows  the  prediction  of  the  outcome  of 
various  experiaental  manipulations.  One  set  of  manipulations  that  can  be  made 
changes  the  intrinsic  frequency  in  part  of  the  cord;  this  is  done  by 
partitioning  the  bath  in  which  the  spinal  cord  is  held,  and  using  the 
dependence  of  the  frequency  on  the  concentration  of  the  excitatory  amino  acid 
that  is  used  to  evoke  the  fictive  swimming  behavior.  The  aatheaatics 
predicts,  for  example,  that  if  there  is  a  local  aaxiaum  in  frequency  that  is 
sufficiently  different  froa  the  other  neighboring  intrinsic  frequencies,  then 
waves  can  propagate  outwards  from  that  point;  thus,  a  stretch  of  cord  can 
"swia  backwards".  One  way  in  which  such  a  local  aaxiaua  is  created  is  if  the 
natural  frequencies  at  a  given  concentration  have  a  rostral  to  caudal  gradient 
(which  is  frequently  seen),  and  if  the  frequency  of  the  caudal  half  is 
artificially  changed.  An  experiment  involving  raising  the  frequency  of  the 
caudal  half  has  been  done  by  K.  Sigvardt  [11],  and  the  backward  swimming  has 
been  observed.  Experiments  involving  directly  raising  the  frequency  in  a  very 
saall  region  so  as  to  create  a  local  aaxiaua  in  frequency  have  not  yet  been 
done. 

Another  experimental  manipulation  for  which  there  is  at  least  some 

preliminary  data  involves  aechanical  forcing  of  the  isolated  spinal  cord. 

/ 

(See  Section  2.)  Numerical  simulations  of  this  forcing  produced  unexpected 


results:  We  found  that  if  H+  and  H  are  multiples  of  one  another  (which 
implies  that  and  rR  are  of  different  signs),  then  phaselocking  was 

difficult  to  achieve  with  forcing  at  one  of  the  ends.  For  example,  if  the 
oscillators  have  uni fora  frequency,  and  the  unforced  waves  travel  in  a  caudal 
direction  with  a  boundary  layer  at  the  caudal  end,  then  forcing  at  the  caudal 
end  (by  any  frequency  different  froa  that  of  the  non-forced  phase locked 
enseable)  appears  to  lead  to  loss  of  locking,  especially  if  the  ascending 
coupling  is  weaker  than  the  descending  coupling.  If  H*  and  H  are  such 
that  r ^  and  yR  are  of  the  same  sign  (e.g. ,  negative,  which  corresponds  to 
waves  in  the  caudal  direction  as  above) ,  then  the  forcing  at  the  caudal  end 
does  not  appear  to  prevent  entrainment,  and  significant  differences  in 
frequency  between  the  unforced  phase locked  ensemble  frequency  and  the  forcing 
frequency  can  be  sustained  without  the  loss  of  locking.  It  can  be  shown 
analytically  [35]  that  if  and  yR  are  both  negative,  then  forcing  at  the 
caudal  end  at  higher  than  the  unforced  frequency  of  the  locked  chain  causes 
the  phase  lags  to  tend  toward  or  beyond  zero,  and  forcing  at  the  rostral  end 
at  higher  frequencies  causes  the  lags  to  become  larger  in  absolute  value. 

(The  opposite  is  true  for  and  fR  >  0.)  This  experiment  has  been  done  on 
the  isolated  lamprey  cord,  and  locking  has  been  found  for  a  significant  range 
of  forcing  frequencies,  with  forcing  at  either  end  [16].  Furthermore,  the 
lags  behaved  as  described  above  for  yR  <  0.  There  was  a  much  larger 
effect  for  caudal  forcing,  which  is  consistent  with  having  a  stronger 
descending  coupling. 

t 

One  way  to  achieve  the  kind  of  anisotropy  in  which  and  yR  have  the 

same  sign  is  to  have  differences  in  timing  between  the  coupling  signals  of  H+ 


T 


T 


and  those  of  H  [31];  such  differences  have  been  found  in  the  lanprey  cord, 
along  with  differences  in  strength,  for  the  part  of  the  coupling  due  to 
excitatory  synapses  [14].  (The  ascending  excitatory  coupling  appears  to  be 
auch  weaker'. )  It  is  also  possible  that  crossed  inhibitory  coupling  plays  a 
larger  role  in  either  the  ascending  or  descending  coupling.  Since  forced 
chains  of  oscillators  with  such  weaker  ascending  coupling  do  not  easily 
phaselock  to  a  range  of  forcing  frequencies  at  the  caudal  end,  crossed 
inhibitory  coupling  nay  be  important  in  the  ascending  coupling.  More  work  on 
forced  chains  of  oscillators  is  in  progress. 

The  anisotropy  suggested  by  the  forcing  experiaents  have  several 
inplications.  One  is  that,  if  H+  and  H  are  thought  of  as  separate  parts 
of  the  coordinating  systea,  each  alone  would  produce  waves  in  the  caudal 
direction.  Another  implication  is  that  backward  swiaaing  in  the  laaprey 
requires  aore  than  a  snail  perturbation  in  the  bias  of  the  coupling.  (See  the 
previous  section.) 

D.  Averaging  and  circuit  design 

In  Section  4F,  we  saw  that  oscillators  coupled  in  both  directions  by 
means  of  signals  localized  within  the  oscillator  cycle  nay  lose  rhythmicity 
and  go  to  a  stable  steady  state,  and  that  this  can  be  avoided  by  having 
coupling  signals  that  are  aore  dispersed  within  the  cycle.  Numerical  work 
shows  that  even  if  the  nuaber  of  coupling  signals  per  cycle  is  small  (e.g.,  3 
or  4),  oscillator  death  can  be  prevented  [28].  A  neuronal  oscillator  with 
many  coupling  signals  can  easily  be  constructed  if  the  oscillator  is 
composite,  i.e. ,  composed  of  a  set  of  neurons  (each  of  which  aay  or  aay  not  be 
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competent  to  oscillate)  and  which  can  burst  at  different  tines  in  the  cycle. 
The  leech  segmental  oscillator  has  this  property;  identified  cells  in  the 
oscillator  have  bursts  with  midpoints,  relative  to  some  arbitrary  zero,  of  0, 
90,  180,  and  270  degrees  respectively  [4];  more  recently  identified  cells 
burst  with  midpoints  near  45  and  240  degrees  [36] . 

As  mentioned  in  Section  2,  the  lamprey  spinal  cord  is  not  segmented; 
there  is  no  anatomical  evidence  of  "unit  segmental  oscillators"  as  there  is  in 
the  leech.  However,  the  above  notion  of  composite  oscillators  (having  cells 
of  many  types)  can  be  generalized  to  be  compatible  with  an  unsegmented  cord. 
That  is,  a  description  of  the  cord  as  a  chain  of  discrete  composite 
oscillators  assumes  that  the  coupling  among  the  cells  of  a  particular  unit 
oscillator  is  very  strong  compared  to  the  coupling  between  units.  If  this 
restriction  is  relaxed,  we  have  a  set-up  in  which  small  portions  are  competent 
to  oscillate,  and  there  is  no  clear  distinction  between  intra-oscillator  and 
inter^oscillator  coupling. 

It  is  therefore  of  interest  to  understand  mathematically  what  indeed  does 
happen  when  this  restriction  is  relaxed.  It  is  possible  to  construct 
composite  oscillators,  e.g.  arranged  in  a  ring,  with  phase  differences  within 
the  components  of  the  composite  oscillator  [37].  Preliminary  work  has  shown 
that  when  such  composites  are  coupled,  with  each  element  coupled  to  its 
neighboring  analogues  with  strengths  comparable  to  the  intra-oscillator 
coupling,  the  entire  ensemble  behaves  like  the  chains  of  discrete  oscillators 
that  we  have  been  discussing.  It  remains  to  show  that  the  this  is  also  true 
for  more  generally  constructed  composite  oscillators. 


E.  Develggggntal  iaplicationa 


The  above  mathematics  has  implications  for  the  embryological  development 
of  motor  behavior  in  fish.  Early  coordinated  motion  called  ”C-coilsn  in  fish 
(and  in  all  vertebrates)  involves  simultaneous  contractions  of  muscles  on  one 
side  of  the  spinal  cord.  Somewhat  later  in  development,  fish  embryos  display 
travelling  waves  (and  S-waves,  which  involve  antiphase  coupling  of  units  on 
opposite  sides  of  the  cord)  (38].  According  to  the  mathematics,  C-coils  are 
to  be  expected  if  the  neural  or  muscular  units  are  coordinated 
electrotonically.  Evidence  for  such  early  electrotonic  coupling  has  been 
found  (39,40].  For  neural  units  coupled  "synaptically",  the  expected  behavior 
is  travelling  waves.  Thus  the  mathematics  suggests  that  the  change  from  the 
C-coils  to  travelling  waves  corresponds  to  the  development  of  appropriate 
functional  synapses.  Furthermore,  since  co-functioning  of  presynaptic  and 
post-synaptic  cells  is  thought  to  be  useful  in  the  strengthening  of  synapses, 
the  initial  electrotonic  coordination  nay  be  an  important  prerequisite  for  the 
later  organization. 

F.  Future  directions 

The  work  described  above  concerns  the  effects  of  local  coupling  within  a 
spinal  cord  isolated  from  effects  of  brain  stem  control  and  sensory  feedback. 
It  was  shown  that  such  organization  is  capable  of  generating  the  most  basic 
elements  of  the  stereotypic  motion,  and  of  performing  certain  kinds  of 
regulation  (maintaining  phase  lags  while  increasing  the  swimming  speed, 
buffering  somewhat  against  the  effects  on  phase  lags  of  differences  in  local 


frequency).  However,  the  phase  lags  in  the  isolated  cord  are  not  constant 
along  the  cord  (they  average  about  lfl!  per  cycle  but  can  vary  froa  below 
0.05X  to  above  2X  per  cycle  in  different  portions  of  the  same  cord.) 

Indeed,  the  mathematics  says  that  unless  the  local  frequencies  and  coupling 
strengths  are  very  carefully  regulated,  such  differences  in  phase  lags  are  to 
be  expected.  In  addition,  it  has  not  yet  been  fully  investigated  how 
efficiently  the  very  reduced  system  can  respond  to  the  kinds  of  perturbations 
that  would  occur  in  vivo. 

It  is  unlikely  that  the  regulation  of  phase  lags  in  the  freely  behaving 
animal  is  performed  at  the  level  of  the  brain.  There  are  two  more  local 
levels  that  have  not  been  taken  into  consideration  in  the  above  analysis.  The 
first  involves  again  just  the  isolated  cord;  in  addition  to  the  multiple,  but 
local,  coupling,  there  is  also  believed  to  be  more  long-range  coupling  [17]. 

We  conjecture  that  the  latter  coupling  may  play  a  role  in  the  correction  of 
inappropriate  phase  lags.  That  is.  if  frequency  differences  are  too  large  to 
be  buffered  by  the  multiple  coupling,  or  there  is  a  transient  perturbation 
that  disturbs  the  phase  lags,  long  distance  coupling  can  measure  and  correct  a 
"mismatch"  between  the  actual  and  desired  lags,  especially  if  the  desired  lag 
between  the  two  points  is  approximately  180°;  this  type  of  correction  is 
difficult  to  perform  by  local  interactions,  since  the  phase  difference  between 
successive  oscillators  is  so  small.  The  role  of  this  long  range  coupling 
during  fictive  swimming  may  perhaps  be  assessed  experimentally,  since  the 
fibers  for  long-range  coordination  lie  in  a  different  tract  than  a  portion  of 
the  short  range  connections,  and  hence  can  be  lesioned  without  totally 
eliminating  short-range  coupling. 


V 


The  other  major  source  of  local  regulation  is  through  proprioceptive 
feedback.  Though  there  is  an  extensive  literature  on  the  hydrodynamics  of 
fish  swimming  [41-44],  until  recently,  there  has  been  relatively  little  work 
on  the  mechanics  of  swimming  and  its  relationship  to  the  electrical  activity 
[44-46] .  The  mechanics  is  linked  to  the  electrical  activity  of  the  spinal 
cord  through  mechanoreceptors  in  the  spinal  cord  [46] .  We  conjecture  that  the 
resultant  feedback  smooths  out  the  irregularities  in  the  phase  lags  seen  in 
the  isolated  spinal  cord  preparation. 
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FIGURE  CAPTIONS 


1.  Schematic  drawings  of  the  solutions  to  (6)  with  H+  =  H  =  A  cos  r  + 

B  sin  r,  A  <  0. 

a)  /J(x)  a  -  p  <  0 

b)  p(x)  =  -  p  +  p^(x)  where  p^(x)  is  a  small  positive  bump  function 
with  a  maximum  at  Xq.  The  added  frequency  has  an  effect  on  phase 
differences  only  on  the  support  of  p^. 

c)  p(x )  =  -  p  +  p2(x),  where  PgCx)  is  a  sufficiently  large  bump 
function.  For  x  <  Xq  and  away  from  the  boundary,  r(x)  >  0,  so  the 
wave  goes  backward. 

The  details  of  the  solutions  depend  on  the  signs  of  f"  and  Q. 

2.  A  pulse  function,  centered  at  9  =  0. 


